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are  determined  from  the  Volterra  integral  equations  given  by  the  strain 
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However,  known  constant  strain-rate  response  does  not  determine  the  kernel' 
Variation  of  the  kernel  within  a  given  qualitative  shape  can  leadifi 
different  shapes  of  constant  strain-rate  response,  so  that  both  constant 
stress  and  constant  strain-rate  tests  may  be  necessary  to  deduce  the  optlmt 
single  integral  approximation,  in  preference  to  multi-step  stress  tests. 

Phase  II:>  A  frame  indifferent  differential  operator  law  relating  stress, 
stress-rate,  strain,  and  strain-rate  is  constructed  to  describe  the 
qualitative  features  of  both  constant  load  and  constant  displacement  rate 
response  in  uni-axial  stress  experiments.  No  lower  order  differential 
relation  can  describe  both  responses,  but  the  two  responses  are  not 
sufficient  to  determine  the  response  coefficients  of  the  relation.  The 
jump  relation  is  determined  for  a  stress  discontinuity  applied  In  a  general 
configuration.  A  simple  initially  isotropic  model  is  proposed  to 
investigate  the  effects  of  loading  history  and  the  anisotropy  induced  in 
a  configuration  following  a  loading-unloading  cycle..  Both  uni-axial  stress 
and  shear  cycles  are  treated,  followed  by  uni-axial  4tress  loading  In 
different  directions.  The  respective  deformation  histories  are  determined 
in  the  small  strain-approximation,  and  demonstrate  a 'significant  influence 
of  the  initial  load  cycle  on  response  from  the  underloaded  configuration. 


PREFACE 


This  report  comprises  two  self-contained  papers 
presented  as  Parts  I  and  II.  Part  I,  presented  as  the 
Fifth  Periodic  Report  in  May  1981,  will  be  published 
shortly  in  Cold  Regions  Science  and  Technology.  It  is 
an  investigation  of  viscoelastic  solid  relations  of 
differential  type  which  can  be  determined  by  constant 
stress  and  constant  strain-rate  responses  of  ice. 

Part  II  is  an  investigation  of  various  classes  of 
viscoelastic  integral  operator  relations  to  see  if 
improved  and/or  more  tractable  models  of  ice  response  can 

be  constructed.  Illustrations  and  conclusions  particular  i 

i 

to  ice  response  are  given.  This  paper  will  be  submitted  1 

to  Mechanics  of  Materials.  Alternative  representations 
to  improve  the  correlation  with  ice  data  have  been 
constructed  in  the  current  period  and  will  be  described 
in  the  next  periodic  report,  February  1982.  An  associated 
paper  will  be  submitted  to  Cold  Regions  Science  and  Technology. 
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Viscoelastic  Solid  koi.it  ions  For  The 
Deformation  Of  Tee 

U.  Spring  and  L.  W.  Mori and 

School  of  Mathematics  and  Physics 
University  of  Fast  Anglia,  Norwich  (U.K.) 

Abstract 

A  frame  indifferent  differential  operator  law  relating 
stress,  stress-rate,  strain,  and  strain-rate  is  constructed 
to  describe  the  qualitative  features  of  both  constant  load 
and  constant  displacement  rate  response  in  uni-axial  stress 
experiments.  No  lower  order  diiferential  relation  can  describe 
both  responses,  but  the  two  rest-  'uses  are  not  sufficient  to 
determine  the  response  coefficients  if  the  relation.  The 
jump  relation  is  determined  for  .1  stress  discontinuity  applied 
in  a  general  configuration.  A  simple  initially  isotropic 
model  is  proposed  to  investigate  the  effects  of  loading  history 
and  the  anisotropy  induced  in  a  configuration  following  a 
loading-unloading  cycle.  Toth  uni-axial  stress  and  shear  cycles 
are  treated,  followed  by  uni-axial  stress  loading  in  different 
directions.  The  respective  deformation  histories  are  determined, 
in  the  small  strain-approximation,  and  demonstrate  a  significant 
influence  of  the  initial  load  cycle  on  response  from  the  un¬ 
loaded  configuration. 


Introduction 


Engineering  application  in  -ninly  concerned  with  the 

mechanical  response  of  polycry s*  ..  1  i  ine  ice  to  maintained 

stress  up  to  strains  of  a  few  pot  cent,  which  includes  the 

instantaneous  elastic  strain,  primary  decelerating  creep, 

secondary  or  approximately  steady  creep,  and  some  part  of 

the  accelerating  tertiary  creep.  i'i  iure  1  shows  the  typical 

strain  and  strain-rate  responses  to  a  constant  uni-axial 

compressive  nominal  stress  o  (force  per  unit  initial 

cross-section )  at  constant  temperature  where  the  longitudinal 

engineering  strain  e  is  the  contraction  per  unit  initial 

length.  Constant  a  corresponds  to  constant  axial  load, 

which  is  a  common  test  situation  (Mellor  1980).  Similarly, 

constant  e  corresponds  to  constant  end  displacement  rate, 

a  common  alternative  test  configuration,  and  the  associated 

stress  response  at  constant  temperature  is  shown  in  Fig.  2 

as  a  function  of  strain  e  =  ot.  Since  the  Young's  modulus 

and  bulk  modulus  of  ice  are  of  cider  lO^°Nm  ^  (Sinha  1978a, 

Michel  1978,  Mellor  1980),  a  moderate  stress  of  lO^Nm  ^ 

-4 

induces  an  elastic  strain  of  magnitude  10  compared  with 
creep  strains  of  magnitude  10  “  common  in  application,  but 
elastic  strains  are  retained  in  the  present  theory.  The 
initial  elastic  strain  is  ee  f  ' )  c/f,  ,  linear  in  stress, 
where  EQ  is  the  Young’s  modulus  at  zero  stress.  The  strain 
jump  induced  by  a  stress  jump  applied  at  an  arbitrary  stress 
level  will  be  derived  as  a  limit  of  smooth  loading. 

In  Fig.  1(b),  em(a),  t  ( i ) ,  are  the  stress  dependent 
minimum  strain-rate  and  time  to  minimum  strain-rate,  and 
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em(o)  in  Fig.  1(a)  is  the  corresponding  strain  at  the  first 

inflexion  point  of  the  creep  curve.  The  time  t  (o) 

in 

increases  significantly  as  decreases,  and,  more  generally, 

so  do  the  times  to  reach  any  r  >in  on  the  creep  curve;  that 

is,  the  rate  process  decrease:  a:  a  decreases.  Thus 

•  — 

em(o)  is  an  increasing  function  .  >t  ',  and  we  suppose  also  thal 

the  initial  strain-rate  e  (o'  increases  with  o.  The  I 

1 

rate  of  decrease  of  lengtli  per  unit  current  length,  or 
compression  strain-rate  with  respect,  to  the  current 
configuration,  is  given  by 

r  =  (l-e)_ie  ,  e  1  ,  (1) 

and  the  limit  e  1  corresponds  to  the  element  length 

shrinking  to  zero.  In  glaciology  it  is  commonly  assumed 

that  there  is  a  long-time  steady  creep  response  r m(o) ,  where  < 

is  the  Cauchy  stress,  which  is  an  increasing  function  and  zero 
at  a  =  0.  With  the  later  o  ,  o ,  e  relation  it  then  follows 
that  .  . 

e  ■*  e  1 ,  o  -*■  0,  r  •>  0  as  t  »  ,  (2) 

but  this  later  part  of  the  tertiary  response  with  strains 
of  order  unity  will  not  be  reached  in  many  engineering 
applications,  particularly  at  low  stress.  We  will  focus  on  the 
small  strain  response  shown  by  solid  lines  in  Fig.  1  and 

Fig.  2. 

At  constant  strain-rate  e  ,  limited  by  e  <  1,  the 
maximum  stress  is  oM(e),  and  the  strain  and  time  at  maximum 
stress  are  eM(e),  tM(e)  ,  where  eV]  =  et^.  As  the  strain 
increases  beyond  eM  the  stress  e 


relaxes,  but  must 


subsequently  increase  a..  <■  -  1. 

that  both  em  and  ev  arc  upp:  a 
range  of  stress  and  strain-rate  t 
there  are  indications  that  the  m 
the  constant  stress  >'  recuiiiM 
strain-rate  e  =  c;  '  hat  ; 

lTl 

em[o  (e)  ]  •-  e  and  i\ 

The  non-monotonic  straia-ra- 
stress  and  non-mono ton ic  stress  ; 
strain-rate,  shown  in  I-igs  i  .u, 
viscoelastic  differential  opor  if. 
(Morland  and  Spring  1981).  The  f 
the  strain-rate  with  respect  to  *. 
the  Cauchy  stress  (force  per  an  it 
that  the  minimal  form  to  describe 
requires  strain-rate,  st.rain-acc< 
stress-rate.  There  is  no  depend* 
configuration,  and  every  con fi gut 
Thus,  the  model  can  neither  doset 
which  may  arise  in  the  ice  format 
later  configurations  induced  by  t 
there  is  no  elastic  (jump)  stiai: 
The  strain-acceleration  term  is  : 
response  to  constant  stress,  and 
necessary  for  the  non-uniform  so 
strain-rate . 

We  now  present  an  analogous 


'■lellor  (1980)  suggests 
■xi mutely  o*ol  over  a  wide 
•espectivcly ,  and  also  that 

i x i m am  stress  o . .  ( e )  is 

M 

‘  '  produce  a  minimum 

i  '  ~  •  (3) 

o  response  at  constant 
o -ponce  at  constant 
2,  have  be*,  n  modelled  by  a 
r  relation  of  fluid  type 
iuid  relation  is  based  on 
he  current  configuration  and 
urrent  area) .  It  is  shown 
-  both  types  of  response 
deration,  stress,  and 
- ico  on  a  reference 
■avion  is  necessarily  isotropic, 
ibe  an  initial  anisotropy 
ion,  nor  anisotropy  of 
lie  loading  history.  Also 
:  when  a  stress  jump  is  applied, 
locossary  for  the  non-viscous 
the  stress-rate  term  is 
ess  response  at  constant 

differential  operator  law  of 


solid  type,  with  explicit  do  i-  -..v  on  the  strain  from 
a  reference  conf  igurut.  Lor. .  -Vm stress  and  stress-rate 
terms  are  necessary  to  ties--!  : -he  response  at  constant 
strain-rale,  but  strain  an  i  i  i  i .  :.-rate  torus  are  then 
sufficient  to  describe  t  to  r  n-  :nif  mi.  strain-rate  response 
to  constant  stress.  That  is,  a  -oust ant  stress,  the  relation 
becomes  a  first  order  d i i Cerent i a  1  equation  for  strain, 
analogous  to  that  for  stra.i  n-itr.s  in  the  fluid  relation,  and 
there  is  an  appropriate  non-un 1  form  strain  and  strain-rate 
solution.  The  significant  mm- 1  i  ueurity  of  the  response 
of  ice  to  stress  even  at.  small  e  rains  means  that  the 
differential  operator  relation  is  highly  nen-iinear.  Linear 
viscoelastic  models  pro post:  1  1  y  Novo  1  ( 1976)  for  low  stress 
applications  offer  well  dove  loped  methods  of  boundary-value 
problem  solution,  but  are  not  ,»r.  appropriate  approximation. 
Sinha  (1973a,  1978b)  and  fold  and  i'inha  (1980)  have  presented 
the  results  of  constant  uni -ax:  :■  stress  tests  for  short  time, 
small  strain,  creep,  and  inferred  an  empirical  relation  for 
the  axial  strain  in  terms  of  stress  and  time  explicitly.  It 
is  not  shown,  however,  how  this  result  can  be  related  to  a 
coordinate  invariant  tensor  law  necessary  for  combined  stress 
loading,  nor  how  the  response  to  non-uniform  stress  history 
is  constructed  in  the  absence  of  linear  superposition. 

While  the  forms  of  the  uni-  mini  responses  determine 
the  minimial  set  of  terms  required  in  a  differential  relation, 
they  cannot  determine  the  shape  of  t:he  tensor  relation  which 
requires  the  response  to  tri-.-xi  il  stress  cr  combined  axial 


and  shear  stress  (Norland  if'7"). 


■".here  are  an  infinity 


of  tensor  relations  which  1* < •«. i •  i  -«  t  •  •  t  he  same  uni-axial 


response.  In  the  absence  i  I  coni  ;  ned  stress  results  we  adopt  I 
a  limited  shape  cf  tensor  iclnt  i  wi  th  limited  dependence 
of  response  coefficients  on  ;t:  ••  and  strain,  and  examine 
its  implications  for  uni-axial  n::;s.  It  is  related  to 
families  of  creep  curves  for  .It  i  :  «  i  nt  constant  stress,  and 
families  of  stress  curves  f'u  !  •  •  .»•  ent.  constant  strain-rates.  1 

In  addition,  the  response  tn  1  <  ■  t  d  i  :• ■  :-uri  1  oad  i  ng  cycles  is 
investigated,  which  demonstrate  -.  >  nut  a  constant  reduced 
stress  response  contains  no  new  ;  nforn.ation  about  response 
coefficients,  apart  from  prv.-iii  ;  imi  t  functions  at 
o  =  0  when  the  stress  is  ■  1  •! .  1  y  unloaded.  However, 

measured  reduced  stress  re.-ao;..  .  ay  provide  consistency 

tests  for  the  model. 

We  focus  or.  a  mode  l  ware:.  .  i  ••  'tropic  in  the  reference 
configuration,  and  deter:;  i-:c  •  :,!  lcronti.il  equations 

describing  combined  shear  n.o  I  ■•xtensions  in  constant 

plane  stress  for  the  small  dr.r:  approximation.  An  initial 
uni-axial  load-unload  :yclo  i  lysed ,  followed  by  alternative 

uni-axial  loads  in  two  ortiu  :•  >n  .« .  directions.  This  demonstrates 
that  if  the  initial  load  is  ro»  •  .  mod  from  the  unloaded 
configuration,  the  response  is  a  .erent,  which  Is  one  feature 
of  the  strain  dependence  absent  a  fluid  model.  A  distinct 
effect  of  the  strain  dependence  is  the  induced  anisotropy, 
or  material  asymmetry,  in  t  he  a-. .  t  ied  eonf  iquration, 
demonstrated  by  the  different  r  .u.onsos  to  identical  loads 
applied  in  the  two  direct  ions.  shear  load-unload  cycle 
is  also  analysed,  whim  ,  in  the  m  1 1  1  strain  approximation, 
is  accompanied  by  negligible  axis!  -.trains  (a  feature  of 

A 
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linear  theory),  and  so  i  mi  in  «.-■ 
axial  loading.  However,  i  i:.  r 
can  remain  constant  or  decay  :  ;• 
axial  loading,  dependin;  on  t  h< 

Differential  tensor  relation 

To  obtain  non -uni  1  or::*  r *.  nr 

and  non- uni  fern  stress  if 

stress,  a  differential  root  ere  i 

cases  to  at  least  a  first  r  ;  < 

respective  response  v.u  in  in  >  r 

analogy  with  t);e  vi  sec  ... ,  •  ic  :  .  ■ 

this  implies  that  stress  o.*.  .-on 

s  *.  r  a  1 :  i  n  ! 

necessary,  and  that^si. ra  i  i  ■  • 
replacing  i  he  strain-r  .«•  e  ar  :  •  : 

of  a  fluid  rel.it ion  wi  i  ••  * 

We  ado:^  the  usual 

tr  D  ---  '  , 

where  the  strain-rate  w i  >. 
configuration  is  the  syrino*.  i  io  *.  * 
gradient  (and  the  rota!  o.r.  i  ■  « 

F  is  the  deformation  era d  i  <  r.  t  . 

tensors  are 


with  principal  invar iantfi 
=  t.rH  •=  trC  ,  K 
K  , 


asymmetry  in  subsequent 
awn  that  the  shear  strain 
i!  leant ly  in  subsequent 
‘rain-tensor  dependence. 

-  i  t ::  ■  at  c--  rstant  stress, 

!  .•trjin-r.iV.:,  in  uni-axial 
:  1 1:  i  >n  ::vmt  reduce  in  both 
:  :  'erential  equation  for  the 
!  and  md  Sprirr;  1981).  By 
1  .  1 1  f  i  e n  ;  *  1  a  rel  ation 

••  •-!••**  e  .r.u:  t  r  :  :  ut  ior-s  are 
1  -  :  :  •*  )  ■>'.:  -re  necessary, 

.  ,  volt  ;  a’  'n  combination 

;  c .  t  n  r  :  :  r  ,  i  n  . 

•  ilit y  a [ ■; ■no; i ration 

,  (4) 

:  ;  '..'1  to  t!  ■■  current 
:  o‘  the  spatial  velocity 

is  ti.e  skew  part),  and 
.hi  Cauehy-Oreen  strain 

(5) 

'  ■’  -  ;  nr  '  ■  '.{  ( trC )  2  -  trC2} 


(6) 


and  the  principal  str.iin-t.is 


:  i  u.'  s  ,irt' 

I:  =  tr D  :  O,  I,".1!’,  -  ,  .  (7) 

Incompressibility  implies  : .  i !  :i  !  v  *  he  stress  devint.or 

S  =  a  -  ~  !  trc  )  I  ,  (  8 ) 

is  determined  by  the  do  format  t.  hi  story,  where  1  is  the 
Cauchy  stress.  An  associated  •  tie- ind  i  f  ferent  deviatoric 
stress-rate  is 

S(1)  =  S  +  S(D  ~  W)  -  i:  -  his  ,  (9) 

where  a  superposed  dot  dor.-  t  os  mat er ini  time  derivative. 

The  principal  deviatoric  at  t  •  :iv  iriants  are 

JL  --  trS  :  O  ,  d ,  i  ’  dot  h  .  ( 10) 

•  ■  - - 

With  the  further  simp:  if;  -  it  i  that  the  stress  and 

stress-rate  tensors  enter  in  a  i.  u  comb  in  it  ion,  and  that 

the  strain  and  strain-rate  ter..;  :  s  enter  as  separate  terms, 

then  the  frame- indifferent  different  ial  relation  described 

above  takes  the  form 

S  +  |  ( t  rS  1)1  l  1 

^  j 

[F  f  (C)FT  -  i(trc  nil  t  :  j  P  +  ■;  2  [  D2  -  -2  I2  1J.  (11) 

Here  f  is  an  arbitrary  symmetric  tensor  function  of  symmetric 
tensor  argument,  and  the  re  spots..;  ■>  eo,.f  f irients  y ,  i>2  are 


and  their 


functions  of  the  invariants 
rates.  Dependence  on  stress 
not  generally  linear,  and 
strain,  and  their  rates  can  .1;  : 
products.  We  will  present  t  !.• 
configuration  which  is  an 
demonstrate  the  anisotropic  r<: 
reached  by  different  loadin’-:, 
of  (11)  is  given  by  roguirir.  : 
isotropic  function  of  h  ,  air 
dependence  on  D  by  frame-;  a.: 


tad  •<  t  ress—  i  a  *  <>  is  therefore 
1  ed  depeneen  on  stress, 
:s.f  ‘  hr 01.  j !i  invariant 
■  in  ilyr;  is  for  a  reference 
1  .;<io  cor.  f  iguration,  and 
r-ponse  from  configurations 
-he  isotropic  reduction 
(!:<■  F  dependence  to  be  an 
•  cad  .•  imposed  for  the 
i  ;  ’  ferenoo ,  hence 


+  :  [  B 


(12) 


where  the  additional  r  .•••  :  i  o  Lents  .  j  ,  ^  are 

also  functions  of  invar  Laat s . 


While  dependence  on  thu  c • r . .  • ; ,  tensor  u;  absent  if 

~  0,  dependence'  ••  •  :  nr.  can  s till  enter  through 

dependence  of  f,  ;  ,  ,  •  on  •  o  r  iin  invariants  K.  ,  K_. 

Furthermore,  such  depen  !■  .-no.  invariants  K  ,  K,  only 

is  sufficient  to  induce  arris  >t  •■••pio  response  from  subsequent 

configurations.  Considc-c  a  -  t,  i  m 'ferenoo  configuration 

reached  from  the  first  by  a  deformation  F  ,  and  let  F 

denote  the  deformation  from  t : :  •  *  second  reference  configuration 
Then 

F  =  F  F  ,  B  -  F  B  F1  , 

~  "°  '  °  •  (13)  i 

; 

Kx  *  tr(CD0)  ,  K-,  =  \  f  K  j  -  M(C!ioCBo)!  : 


flippy—  - 

-  10  - 

while  the  strain  invariant:;  f  i  on  the  second  reference 
configuration  are 

K :  =  trC,  K.,  o'  K“  -  •  i  "1  .  (14) 

Thus,  unless  D  =  b  1  ,  wla.-i.  ‘  i;e  change  of  configuration 


is  simply 

an  isotropic 

liil.lt  . :  •  i  • 

v>r;pr( 

1  l'r>  S  \ 

on ) 

,  K 

1  '  K2 

cannot  be 

expressed  i n 

O  I.V.J  *.  ! . 

. invar 

l  an 

l 

Kr 

alone 

but  depend  on  the  slr.u 

:i  te:.:;  ; 

c  i.  t 

1  y . 

That  is, 

response 

functions  . , 

;  1  ’  ;  ' 

'!  Kl' 

K? 

OX 

hi  bi 

t 

explicit 

dependence  on 

the 

t.  one. or 

C 

f  r 

cm  a 

second 

reference 

conf igura t ion 

in  goner. 1  1 

,  not  j 

us  t: 

on 

the 

invariants  K^,  1^,  so 

the  second 

con f iqu 

.  r  a  t. 

i  on 

i  s 

not 

isotropic 

when  =  (o 

2  ~  0 • 

■h  anise 

•  tropy 

will 

be 

demonstrated  by  example. 

Jump  relations 

i 

It  remains  to  complement  t1.  11 fferentiai  relation  (12) 

by  jump  conditions  which  dot  i  i  m  i t  he  strain  jump  when  a 
stress  jump  is  applied  in  any  or n !  i cur at  ion  during  loading. 

If  the  jump  relation  is  the  limit  of  a  sequence  of  smooth 
rapid  changes  governed  by  the  a i f re ten t ia 1  relation  (12), 
then  the  response  coefficients  1 ^ 

must  remain  bounded  as  stress-rater,  and  strain-rates  become 
infinite,  which  is  a  restriction  on  the  dependence  on  J2#  Jj* 
and  K^,  I^,  or  I~  ,  T  .  Further,  since  / D  dt  is  unbounded 
as  the  time  interval  approaches  zero,  must  approach  zero 
sufficiently  fast  as  stress  and  strain-rates  become  infinite, 
which  implies  that  a  non-zero  tf>2  must  depend  on  the  invariant 
rates.  Since  uni-axial  response  does  not  seperate  the  and 
4> 2  terms,  andsince  we  subsequently  analyse  models  which 
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exclude  dependence  of  the  response  coefficients  on  the 
invariant  rates,  we  henceforth  adopt  the  simplification 

<p2  =  0  .  (15) 

Suppose  the  deformation  is  changed  instantaneously  at 

time  t  from  F  (t)  to  F  by  a  stress  change  from  S  (t)  to  S. 

^  o  ~ 

Since  such  elastic  strains  are  infinitesimal,  and  assuming 
there  is  no  allied  rigid  rotation  jump,  we  will  adopt  the 
linear  approximations 

W  =  0,  F  +  (1  +  c)  r 

'  ~  ~  '  ~  '°  7  (16) 

•  «.]  •  —  X  • 

D  =  F  F  ,  D  -  FF  =  D  +  c  , 

~o  ~o~o  '  -  o 


(18) 
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where  ip,  <p ^  are  evaluated  in  the  infinite  rate  limit  (if 
dependent  on  invariant  rates)  and  at  F  =  F  ,  but  retain  their 
dependence  on  S.  Given  SQ,  Fq,  the  six  integral  relations  (18) 
hold  for  arbitrary  S  and  determine  the  six  components 
g^j(S,  SQ,  Fq,),  and  hence  h(S,  Sq  ,  Fq)  by  (17)  and  the 

initial  condition  h(S  ,  S  ,  F  )  =  0. 

~  ~o  ~o  ~o 

Now  jjglj  ^  represents  an  elastic  modulus  magnitude, 
so  that  |  |S|  j  .  |  |g  |  |  'v-  j  j  c  (  j  1,  so  consistent  with  the 
infinitesimal  elastic  strain  approximation,  (18)  becomes 

S 

-  «i5ikdskJ>  *  0  •  '19) 

S 

A  solution  of  (19)  is  illustrated  for  uni-axial  stress 
loading. 

Uni-axial  stress 

For  incompressible  material  the  nominal  stress  a  and 
Cauchy  stress  o  are  related  by  (Chadwick  1976) 

o  =  F  a  .  (20) 


stress  they  become  function:'  two  variables  o  and  e  , 
and  explicit  notation  for  t  he  t  <•  h:  rod  dependence  is  useful; 
thus 


i|»  =  ii>  (o  ,e)  ,  t  j  (■’  »t>)  ,  .  ,  ,  (■•  ,  e )  ,  ..  w2  (f  ,e)  .  (27) 

I 

Given  dependence  on  the  i  nvnr  i  ers  detorni  nor.  the  reduced 
functions,  but  the  latter  oam  *  d  i  a t i nauish  dependence  on 
J_  and  J, ,  and  on  K,  and  K  .  .  That  is,  there  are  an 
infinity  of  response  coefficient  sets  compatible  with  the 
same  uni-axial  stress  response,  an  1  the  response  for  two 
independent  stresses  and  two  i nan:  undent  strains  is  required 
to  construct  the  general  ] oad i n  :  re  I  at  ion  for  an  incompressibli 
material.  "or  illustration  .d  : . > •  i r  response  we  interpret 
a  model  dependence  on  ',o,  in  terms  of  j. ,  K ^ .  Clearly 
a  non- symmetric  response  in  .'w::sion  and  tension  would 
require  dependence  on  both  even  and.  odd  invariants. 

Each  principal  component  of  (1.1)  -lives  the  same  relation, 
which,  using  the  restriction  (id  !,  is 

(l-e)3o  +  (1-e)  2[  (l-e)e  -  lie?  1  -  yt^d-eje  x  ue  ,  (28) 


where 


oe  =  u.  (1-e)  f  1  -  ( 1  -o )  3  1  +  -  (i-el('’, 


<•  6j.i0)e  +  Ole") 


(29) 


Only  the  combination  of  occur?;,  so  that 

uni-axiai  stress  response  carrn  >  o  i  st  i  ngu  i  s'i  the  contribution 
and  u' 2 


of  the 


term  ',  i  i>  i  .  1  . 


With  the  restricted 


1  s 


dependence  (27)  the  response  e  ■  \ .  :  icients  ere  necessarily 

bounded  as  rates  become  infinite,  and  the  jump  relation  for 
a  stress  jump  oQ  to  ■  with.  •  responding  strain  jump 

A  _  _ 

e  to  e  +  h(a,o  ,e_)  is  obtained  by  integrating  (28)  with 

O  O  .  O  t  ’  ...  .  .  i 

respect  to  time  over  a  van;  sh .  r.u  time  interval 


a 

J  U  (l-e)3  -  (  -|v  1  i  !  -o !  . 


o 


(30) 


Since  (30)  holds  for  all 


1  !  ::V' 


,  the  integrand 


vanishes  to  give  a  first  a  do?  •.  i  "f  event  ia  l  equation  for 

t\  __ 

the  strain- jump  h  ( a  1> 


(2o(l-e  -h)  t|?  (o  ,  e  +h )  ?  -1:  . 


i  o 


‘!l,l  ' 


A  ^  \ 

=  (1-e  -h)  C'  (o  ,e  +!i)  , 
o  o 


. ;  i  0  fC  :=  0  . 

o  o  o 


Applying  the  infinitesimal  strain  jump  approximation  1 

^  i_  1 

1 h |  <<  1,  \o  — |  <<  1,  and  evaluating  the  response  coefficier* 


at  Fq  ,  that  is  0,  a 

gives  the  simpler  equation 


(even  when  e  -  0) , 
o  o 


3 :  ,  3h  =  (1-e.  1  -  .  1  '  ,<•  )  ,  (32) 

A. 

in  which  h  does  not  appear  in  arguments  of  the  response  j 
coefficients.  While  e  nay  be  neglected  in  comparison  with 
unity  at  small  total  strain,  the  ,  ,  c ^  dependence  on  eQ 

may  be  significant.  The  initial  strain  jump  is 

a  _  _  ■ 

e  (o)  =  h(c,0,0),  a.  2  :t,  (0,0)  ,  1 (0,0)  us  a  -  0.  The 


condition  (32)  is  given  directly  by  (19)  since  h(o,oQ,eo) 

-  -<l-e0)  hn<5-  !0-  ?„>• 
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At  constant  load,  >  o,  reduces  t  an  explicit 

•  •  • 

expression  for  e  in  terms  u  >nd  e;  that  is,  a 

first  order  differential  egu.it.  .  •>..  :  .  .r  o.  din  ce  «  ( t )  is 
monotonic  for  a  >  0  ,  ii;.  In,  the  strain-fate  response 
shown  in  Fig.  lb  can  be  expressed  as  a  function  of  e  instead 

•  j 

of  t,  and  hence  the  family  oi  i  vs  p  ir.se  curves  for  different 

| 

constant  o  determine  a  relat  i  o  ----  f(  -,e)  for 

_  i 

e  s  e  (o),  o  >  0.  By  (281,  1 

6  i 

1.  .  „  ,  ,  1  ■■  _ 

0  =  0:  e  =  F(o,e)  -  - - - - -  51 - ,  e(0)  =  e  (o).  (33) 

3  a  2  e 

j;  j  d-o}  i  2.  '  (1-e)  j 

j 

! 

At  constant  strain-rate  o  -  w,  the  typical  stress-strain  i 

curve  is  shown  in  Fig.  2.  The  ;•  •  ik  stress  r  increases 
with  w,  but  it  was  noted  that  Liu-  corresponding  strain 
eM  is  approximately  constant  i  '  0*01),  so  it  is  reasonable 
to  assume  that  at  each  strain  c  fat  least  in  the  small 
strain  range)  the  stress  increases  with  w;  that  is,  the  curve! 
for  different  w  do  not  intersect  .  Hence  at.  each  value  of 
e  there  is  a  monotonic  ;  -  w  relation,  which  can  therefore 
be  inverted: 

e  =  w  =  constant:  o  =-  ( 1  ( w ,  <:■ )  ,  w  -  W(o,e).  (34)1 

Thus  the  family  of  response  curves  G(w,e)  lead  to  a  family 
of  generalised  Young's  moduli  (  al/iie)  which  can  be  expressed* 

as  a  function  E(o,e).  From  (28),  and  using  the  defenitions  1 

(331  and  (34)  of  F  and  W:  1 


e 


E  ( o  ,  e ) 


(35) 


=  w : 


3  c 
3e 


w 


1/ 


l  i  -  ij 


The  initial  Young's  modulus  (■  •  .  ,<  -  O,  u  -  0)  is 


E 


o 


E  (0,0) 


(36) 


since  F(0,0)  -  0,  which  is  ind.  eondent  of  w.  This  is  a 
consequence  of  assumin'!  rate  in.'’*,  pendent  response  coefficients, 

but  is  a  good  approximation  over  a  practical  range  of  w  with 

0  -  ■>  — 
the  value  EQ  -  9  *  10  Nm  “  iMichoi  1978).  While  E  >>  o, 

E(o  ,e  )  =  0. 

”  _  from  the  uni-axial  stress  tests. 

Given  E(t',o),  W  (c  ,e)  ,  and  •  v  ,e)  ,  /  the  constant 

strain-rate  relation  (35)  dote rn  ;  nos  the  ratio  ,  then 

the  constant  stress  relation  (31)  determines  the  combination 
[(l~e)J  a  -  me] /-p  ,  so  that  a  fu :  t her/ resul t  is  required  to 
determine  both  u  and  v .  In  the  malogous  viscoelastic 
fluid  model  the  equivalent  results  to  (33)  and  (35)  were  not 
completely  independent,  and  the  adopted  model  would  not  be 
compatible  with  arbitrary  !'  and  .  If  at  a  time  t^ 
during  constant  stress  i  loading,  the  stress  is  partially 
unloaded  to  o  (  >  0)  and  then  again  maintained  constant, 
the  strain  undergoes  an  elastic  decrease  to  e(t^+)  governed 
by  the  jump  relation  (32),  and  then  lor  t  -  t^  is  governed 
by  the  differential  equation  (33)  with  a  -»  and  e(t^+) 

as  starting  value.  Given  that  it  ,e)  has  been  determined 
over  the  necessary  range  of  (  ' , o :  by  families  of  loading 


j 
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responses,  then  F(o  ,e)  is  :.r.  ■■•.a  ind  the  response  in 

t  >  is  determined;  that  is,  the  same-  combination 

3  — 

[  (1-e)  o  -  oe]/i p  arises,  and  t  unloading  response  does  not 
separate  w,y' 

However,  in  the  complete  :  >,»-iing  case,  u  =  0, 

t  >  t  :  e  »  F(G,e)  =•  -  r-i-iLLL^l - ,  (37) 

il-e> 

which  determines  F(0,e),  hence  ( t),  e)  /■: .  ( 0,  e)  ,  not 

determined  by  initial  loading  0  (when  o  0) .  We  are 

not  aware  that  details  of  the  response  to  complete  unloading 
have  been  determined.  Throe  pose- hie  general  situations  are 
illustrated  in  Fig.  3;  u)  no  ’■  l  uxation  (e  --  0)  ,  (ii)  partial 
relaxation  (e  <  0,  e  *  et)  •  0)  ,  ,  ;  i  i )  complete  relaxation 

(e  <  0,  e  -*  0)  ,  based  on  the  assumptions  .*>  (0,e)  >  0, 
oj(0,e)  5  0.  It  is  also  feasible  that  creep  continues  (e  >  0)  , 
either  indefinitely  to  a  limit:  or  to  some  maximum  followed  by 
relaxation,  for  example  if  m(o,c)  o  for  all  0  <  e  <  1, 
or  for  some  e  before  chang ing  sign  as  e  increases.  This 
does  not  appear  to  be  a  smooth  limit  of  the  response 
coefficients  for  o  >  0.  Case  (i)  is  equivalent  to 

A 

w(0,e)  =0  at  e  =  e(t^+) ,  or  . (o,e)  so  if  there  is  no 
relaxation  from  any  pre-loaded  conf iguration. 

Idealised  response 

The  reduced  response  coefficients  are  defined  by  the 
functions  F(cr,e)  ,  E(o,e)  reproson!  ing  families  of  responses  | 
illustrated  in  Fig.  1  and  Fig.  ,  together  with  an  extra 
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assumption  to  uncouple  and  .  We  will  now  construct 

a  simple  model  which  has  the  np;  sopr Late  qualitative  features. 

First  introduce  a  reduced  time  *  -  t/'t  (a)  and  suppose 

in 


■j  lus  e 


m 


{.:■),  are  independent 


that  (de/dt)  at  constant 
of  a.  Setting 

em  =  °*01'  =  =  kf(e),  f(em)  =  lf 


(38) 


then  k  is  a  constant  and 

k  =  t  (a)  e  (c)  ,  F  (a  ,e) 
m  m 

tmG)  = 


=  e  =  e  (u)  f  (e)  , 

e”1  (a)  /,nf-1  (e '  )  de 1  , 
m  o 


(39) 


where  the  universal  function  fie!  has  the  same  shape  as 
F(a,e)  for  each  constant  v.  V.'e  take  the  Colbeck  and  Evans 
(1973)  power  law  for  the  minimum  strain-rate: 


e  (o)  -  0*  2 1  o  +  0*  .14  c  +  (,d,T  ,  c  •  lobNm  ,  (40) 

r;  _  -? 

where  the  stress  and  time  units  tre  L0  Nm  “  and  one  year 
7 

a  =  3*15  x  10  s  respectively.  Thus  F(0,e)  -  0,  and  so  there 

is  no  relaxation  after  complete  unloading,  and  case  (i)  of 

Fig.  3  always  applies.  Also  (39;  extends  F  into 

0  s  e  <  e  (0 )  . 
e 

The  main  features  of  f  (e)  shown  in  F’ig.  1  can  be 

represented  by 

f (e)  =  R  +  Be~be  +  Ce“ce  (41) 


where  R  is  the  ratio  of  the  strain-rate  at  "large  e"  (long 
time)  to  the  minimum  strain-rate  (here  independent  of  0). 
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B,b,C,c  are  determined  by  (38t  ,  f’(em)  =  o,  f"  (e L)  -  0, 

where  is  the  strain  at  i u I  1  ex  Lon,  and  f  (O)  =  1  , 

where  f  is  the  ratio  oi  the  initial  strain-rate  to  the 
o 

minimum  strain-rate  (here  in  :<•;>. ••mUmt  of  )  .  Assuming 
R  =  2*5,  suggested  by  Stein.'::  »:u.  ’  ( 1  '>'.>&)  it  if.  i ,  and  making 

the  arbitrary  choices  :  •  -  •  ,  f  then 


R  =  2*5,  B  =  -  5*6,  l.i  -  '  1  !  •  ’ ,  ■  274  -  7  . 

The  reduced  time  t  is  related  to  the  strain  e  by 

’e  -1  .  __  '  fGm  _] 

f  (e')de'  =  em(o)t  =  t  I  f  (e')de',  =  0.00552t 

0  0 


(42)1 


(43) 


for  this  model. 

Figure  4  shows  the  above  f  an  a  function  of  e/e  and  as 

m 

a  function  of  t,  and  the  associated  strain-reduced  time  curve 
F(o,e)  is  now  prescribed. 

It  is  convenient  to  express  the  generalised  Young's 
modulus  E(o,e)  in  the  factored  form 


E (o , e) 


E 


{■ 


(44n 


where  the  factor  E(u,e)  is  o  :•  ,i  m  to  all  strain-rates  w, 
and  E  (0,0)  -  E  .  We  now  exai:'.  ■>  the  stress-strain  response 
at  constant  strain-rate,  ■  -  c.  <  ,  under  three  different 

models  for  E,  all  independent  >d  a: 


(i)  E  =  F.  ,  (ii)  E  *•  !■:  e'i<,iv  +  E,  (iii)  E  =  K  e  800e  ♦  E. 

o  o  1  ,  o  1 


E  =  9  x  lo4  •  l()5Nm  ",  b,  -  20  »  nv'Nm  2. 

O  j. 
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In  the  exponential  decay  case;-.  <  i  i  j  and  (iii),  the  constant 


Ei  (  <<  Eq)  ensures  tii.it 


remains  positive 


at  strains  of  order  o*.y> 

,  area::'  > 

:.  :  has  re  I 

axed  below 

5-2 

10  Nm  ,  so  that  3h/~  • 

■  < '  .  : . 

•  ••  lump  ioI.it 

ion  ( 32) .  At 

constant  e  =  w,  by  ( 3 r- ) 

,  •.  A  4  )  , 

1  t 

|f  =  E(e>{  1  -  '■« 

Lvi 

w  ' 

,1  ,jt  *  • 

i  \ 

/ 

which  has  been  integrals.' 

:  : ;  u : :  i  • 

.  y  ter  '  ;  • 

ibove  F  and 

different  E(e),  yioldi:. 

the  • 

;  -at  ra  i  r. 

rves  shown  in 

Fig.  5  for  four  values  o 

•  :: -.  -in  '.-.1  i  :  1 

decay  of 

case  (ii)  shows  little  d 

:  f  f  ct'e-i 

'  :  < n  t  h"  car 

'  e  r:o, 

in  e  <  0*(-6,  but  the  no 

re  r.i:  . 

■  ■  : y  of  ;  a  . ;  <  • 

1 i i i 1  fails 

to  show  a  required  peak 

at 

*'  m  * 

Both  case  (i)  and  ( 

i  i  ) 

'  v •  requ i re  i 

qualitative 

response,  which  is  relatively 

■  -n  •  i  t  i  ve  to  r 

his  variation 

a 

of  E(e),  so  the  constant  sires. 

:  lactic  n  }•'(■; 

,e)  also 

controls  the  constant  r.train-ra-. 

•  it-sponse.  S 

ince  F(o,e) 

defines  a  family  of  non- 

i n  to  r sec 

t 1  a  i  strain-rate  functions 

of  e  with  ( 3  F/3o )  >  0 

,  there 

:  s  a  unique  curve 

A 

o  =»  o(w,e)  for  each  w 

"•  0  or.  . 

ich  F  -  W, 

and  F  <  W 

for  a  <  o  .  Thus  the  peak  point 

!  -  G  (w 

,em)]  occurs 

m 

on  o  =  o(w,e),  since 

Or./3o) 

">  at  e  -  e 

M,  and 

(3G/3e)  <0  as  a  >  o 

.  The 

'li  ves  o  -  a(w,o)  are  not 

distinguishable  from  the 

SO  lilt  i"! 

r.  a  -  G  !w,e) 

for  E  =  E 

(46 


shown  in  Fig.  5.  In  all  cases  • 
e  -»  0  and  the  initial  stress  ri 
on  the  figure  scales.  Because 
and  (ii)  has  little  influence  o" 


;  1  opt.'  approaches  Eq  as 


case  E  =  E 


is  adopted  for  fut 


<■  aopears  discontinuous 

vatiation  between  curves  (i) 
t  h  e  m<  >d  u  1  u  s  ( o  ,  e )  ,  the 
:  ‘ie r  i  1  lust  r a t ion . 


_  2  ?  — 

ami  (40) 

Since  F(0,e)  =  0  by  (30),,  implying  no  relaxation  after 
unloading,  the  extra  assumption  on  o,  y  must  give 
w(0,e)  =  0.  For  illustration  we  choose 

euj(o,e)  =  a(l-e)"*o  ,  O  i  1 .  (47) 


When  a  =  0,  o  =  0  and  only  strain  invariants  enter  the 
response,  while  for  a  order  unity  the  strain  tensor 

A 

term  we  contributes  equally  with  the  o  term 

to  the  strain-rate  (33).  Now,  by  (33),  (39)  and  (44), 


¥ 


(l-a)o 


Eoem(o)f (e) 


-2(--  1 }  (-  [  1 - — —  ],  (48) 

3e  (;)f(e)  E  (1-e) 

m  o 


both  positive  and  bounded  as  1  •  o  and  e  *  0  . 

elastic  strain,  setting  eo  -  c  in  (32),  is  ee(a) 
To  obtain  response  coeffioien's  ^ u  as 

of  the  stress  and  strain  invar  i  a:.  >•  s  for  multi-axial 
load  illustrations,  we  assume  dependence  is  on  J2 
only,  and  use  the  small  strain  approximations  e  << 


The  initial 

-  ?/*„• 
functions 

and  shear 

and 

1 .  From 


(24)  and  (26), 


e 


* 


(49) 


then  (40)  becomes 


e_  (a)  =  o(0*21  +  0  •  4  2  J 4  ()*4->‘3d0; 
m  i  j. 


~  od(J2) 


(50) 


i 
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Choosing  =  0  in  the  definition  (29)  of  w  and  neglecting 
e  compared  to  unity,  the  model  (47),  (48)  gives 


w. 


Kr3 


'  Eo* 


1- 


fi(J2)fL( 


(51) 


^1  =  |eq0  -  4(J2/3)!*, 

where  f(e)  is  defined  by  (41),  (42),  and  the  simplified 

relation  (12),  with  (15),  is 

S  +  4i[S(1)  -  |tr(SD)  11  =  +  -o L [ B  -  ^Kxl]  .  (52) 

Loading  cycles  and  induced  anisotropy 

To  demonstrate  the  anisotropy  induced  in  configurations 
-loading 

reached  by  a  stress^-un Loading  cycle  from  an  initial  isotropic 
configuration,  we  will  determine:-  md  compare  strain  responses 
predicted  by  the  relations  (3!),  for  i  =  0  and  a  i-  0 

when  the  initial  cycle  is  on  uni-  ixiil  stress  and  when  the 
initial  cycle  is  on  shear  stress.  r.i  both  cases  the  initial 
constant  stress  is  applied  over  a  time  interval  (0,t^), 
where  t  =  t^  is  in  the  tertiary  creep  range,  and  unloaded 
completely  at  t  =  t^  for  a  time  interval  (t^,t2) •  Since 
the  model  (51),  (52)  is  non-relaxing,  the  strain  remains 
constant  over  (t^,t2).  At  time  t,  ,  two  alternative  reloading 
situations  are  considered  to  compare  the  induced  directional 
properties  of  the  configuration  at  time  t2;  namely,  the  same 
constant  uni-axial  stress  for  t  •  t2  applied  (separately) 
in  two  orthogonal  directions.  Table  1  shows  these  four  loading 
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histories,  further  illustrated  by  the  section  (a)  of  Figs  6-8. 
Non-symmetry  of  the  nominal  stress  implies  ¥  a2i 

in  shear,  and  zero  normal  traction  on  the  sheared  faces 

=  constant  requires  one  normal  nominal  stress  non-zero,  here 


21  "  V  °  1 2  '  °  1 1 

other  o  ^  j  -  0 


/ 


11 


-  0. 


o 


22 


t  >  t 


2  : 


other 


0 


other 


i] 


=  0 


Table  1.  Lcading-unloudina-re loading  cycles 


The  above  examples  invoLve  i  n  i  y  plane  stress  in  OX^X^* 
and  the  principal  stretch  along  ox^  is  given  by  the 
incompressibility  condition.  In  addition,  only  the  (1,1),  (1,2) 
(2,2)  components  of  (52)  are  independent  equations,  so  only  \ 
0X^2  components  of  the  various  tensors  are  displayed.  For 
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combined  stretches  A^  along  OX^,  ^aloivj  ox.,  and  shear 
tan-^- (-^-)  between  OX^  and  OX2,  the  general  forms  are 


h  V 

1 

0  1 1  0  i  2 

.  0  XT  . 

,  F33  ('i‘2-  '  " 

0 21  °22  , 

(53) 


3  l 

°11 

0 12  ]  = 

Vlln  '21 

'  1  °  1 2  ^  X  e  2  2 

1 

°21 

°22  J 

t  X2C21 

A2°22 

using  (20),  where  A.,/^  are  principal  stretches  only  if  the 
shear  y  =  0.  By  symmetry  of  the  Cauchy  stress 


1_ 

\ 

1 


/ 


°21  =  °'  Y  -  0 


\  ^2  - 
X1  21 


>  1 


n 


(54) 


The  latter  situation  is  adopted  for  the  shear  cycle,  so  that 
no  normal  load  is  applied  in  eon  iunot.ion  with  the  nominal  shear 
stress  02p  on  the  loaded  plane--  -  constant  (X2  =  constant) 
Surface  traction  is  given  by  (Chadwick  1976) 

t,  .da  =  oTn  da  =  aT  NdA  ,  (55) 

where  n  da  and  NdA  are  the  current  and  initial  elements 
of  vector  area.  On  a  sheared  face  X^  =  constant,  I 
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2 1  - 1/« 

N  =  (1,0) ,  n  =  1  +  ( 1 , 

■  I  »; 

2 

aTN.n  «  (an  -  ^  *21  «  • 


so  for  y  °21  ^  zero  nom.il  t  r.i;tion  requires  that  either 

°11  °r  °2 2  ;'-s  non-zero»  anci  avi  1  ting  (54)  , 

°11  =  °21  '  1 2 2  "  °  •  (57) 

It  is  found  that  neither  \  nor  '  can  remain  constant  during 
the  above  shear  loading  governed  by  the  relation  (52) . 


Corresponding  to  (53),  noting  D  +  W  =  F  F  t 


2  1 

XJ+Y  y\2 


X2  x* 


f-'I1  l  2)~11  B33  =  (X1X2)-2 

i  I  r (58) 


- 1  : '  -1 

• .  I  i\ ; 


F33  "  X1X2 


?  “  1 


2X1o11-X2022  +  2y°21  3  X  2 2 1 


3X2°21 


2X2°22~X1°1  r>02lj 


'  S33  =  ~(S11+S22) '  (59) 


D  +  W  = 


>IlJi  b1;-(uy2>'l!i 


x-h2 


,  (I)  *  W) 


(x1x2)  * 

~iV 


2  2  -2-2  2 
K1  *  X1  +  X2  +  X1  X  2  4  Y  '  K2 


2  ?  -2  -2  2-2  -2 
12  1  2  Y  Ai  ' 


1  2 


For  small  total  strain,  )  e2  -  1-Xj,  e 2  -  l-x^,  Y^ 


3  3  3 

neglecting  e.^,  e2,  y  , 


*1  "  K2  »  3  +  4(e2  +  e2  +  e^)  +  y2  , 


so  that  dependence  on  and  on  cannot  be  distinguished 
to  this  order  of  magnitude.  Further,  the  shear  condition  (57) 
to  lead  order  is 


The  small  strain  approximation  is  now  applied  to  each  of  the 
four  loading  histories  shown  in  Table  1. 

During  the  initial  constant  uni-axial  stress  =  -  oQ  , 


X2 

-  X  3 

rtf* 

1 

II 

o 

II 

>- 

e2 

1  > 

R11 

=  1  -  Ze.^, 

K1 

-  3  = 

3e2  , 

S11  -  - 

2— 

T'o  ' 

J2  - 

1-2 

3Jo 

'  Dn  ei' 

W  = 

0, 

tr  (S  D) 

=  °o«l' 

(  c;  p ) 

1 r  !r ;  1 1 

-  ~o 

0^1  ' 

®n  =  Kev  and 

(52)  becomes 

e1(3<>1  +  4 voQ)  =  2aQ  -  6^^,  e}  (0)  =  aQ/E0  «  lo'2,  (65) 

which  is  the  small  strain  approximation  of  (33) .  With  the 
response  coefficients  (51),  and  2 oQ  -  6u>^e^,  each 

have  a  factor  (1-a),  so  the  differential  equation  (65)  is 
independent  of  a  ,  and  integrating  over  (0,t^)  gives  the 
same  e^(t),  e2(t)  =  -*ie^(t)  for  all  a  .  Unloading  at 
t  *  t^  gives  a  negligible  elastic  strain  jump 


28 


[e^]  =  -  00/£q'  [e23  =  *ic0/i:o  '  fallowed  by  constant  ei*e2 

over  (t.,t~).  Reloading  to  u ,  ,  -  -  c  for  t  >  t_ 

L  Z  Z 

recovers  the  elastic  strain  !  1  •=  -'Q/^0>  I  62]  =  -!50o/Eo' 

and  continues  the  response  as  if  unloading  had  not  taken  place, 

since  the  same  differential  equation  (65)  holds  with  initial 

condition  e^(t2+)  =  e  (t^-)  ,  •  «'*  .  Figure  6(b)  shows  the 

-  5-2 

strain  histories  e^(t),  e9(t)  when  o  =  10  Nm  and 

t^  =  200  hours  compared  with  t  (  To)  =  119*4  hours.  Since  the 

initial  loading  passed  the  minimum  strain-rate  time  t  ,  the 

m 

response  from  the  unloaded  state  does  not  show  a  minimum  strain 
rate,  nor,  of  course,  any  of  the  initial  leading  features. 

This  is  a  crucial  feature  of  a  viscoelastic  solid  dependence 
on  a  reference  configuration:  identical  loading  tests  on 
samples  with  different  pre-loading  histories  from  the  reference 
configuration  exhibit  different  responses. 

A  distinct  feature  is  that  of  the  induced  anisotropy, 
or  material  asymmetry,  in  the  new  unloaded  configuration,  even 
when  the  material  is  isotropic  in  the  reference  configuration. 
Consider  the  alternative  uni-axial  stress  loading  ^  22  ~  ~  °o 
for  t  >  t.,.  There  is  an  elastic  strain  jump 

Ce2^  “  <W  ^el^  =  aL  t  ~  t2 '  so  that 

e2(t2  +  )  =  e2(tl-)  +  §VEn'  “  _  ^l  (66) 

since  an  <<  j E  e_ | .  In  t  >  t0, 
o  1  o  2  1  2 


e!  “  e*  =  -  *5(e2  +  ^e*)  ,  Kl  -  3  -  3e2  +  Je*2  , 
e2(34>1  +  4iJ;oo)  =  2oQ  -  . 


(67) 
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While  (67)  ^  is  the  same  differential  equation  for  e2(t)  in 
t  >  tj  as  the  above  reloading  equation  (65)  for  e-^(t),  the 
equations  for  the  corresponding  axial  strain  changes 
y2  ~  ®2  4  ^el  anc^  =  ei  "  c*  are  not  the  same,  and  the 

invariant  K^-3  entering  the  coefficients  is  a  different 
function  of  y2  and  y ^  respectively  in  the  two  cases. 

By  (65) 

Kl"3  =  3(y2  +  Sv^e*  +  e*2),  V x  ( 3 1  x  +  4,'<tq)  =  2oq  -  6o)1y1  -  6a>1e£ ,  (68) 
with  y^(t2)  =  (->t  an<^  ky  (67), 

Kj^-3  =  3  (y2  -  y2e*  +  e*2),  y 2  ( 3 1  +  4:f  iq)  =  2oo  -  6u1y2  +  30)^*,  (69) 

with  y2  (t2)  =  0.  K^-3  =  3e£2  at  t  -  t2  in  both  cases 
but  has  different  dependence  y^  and  y2,  so  the  differential 
equations  for  y ^  and  y2  are  not  the  same,  even  when 
0)^  =  0  (a  =  0).  Thus,  y2  /  y ,  in  t  >  t2,  and  the 
material  is  exhibiting  a  different  response  to  loads  applied 
from  the  unloaded  configuration  in  the  0X2  and  OX^ 
directions.  Figure  7(b),  (c)  show  the  results  for  a  =  0 

and  a  =  0*5,  comparing  o2(t)  in  t  >  t2  for  the  a2 
loading  to  e-^(t)  =  e^(t)  for  n  ^  loading  relative  to 
the  same  starting  value  at  t., ,  that  is,  comparing  y2  and 
y^.  The  difference  y2~y^  reflects  the  degree  of  induced 
anisotropy.  The  effects  of  changing  a  ,  hence  ,  are 

not  as  marked  as  the  K^-3  change  in  the  arguments  of  j 

and  ^  ,  because  of  the  strong Lv  non-linear  function  f. 

In  linear  viscoelasticity  all  coefficients  are  constant,  but 


ice  exhibits  strongly  non-line. n  response  even  at  small 


aQ  <<  I EQe | ,  and  the 


strains.  For  the  model  (51),  0!  !tj_cl  since 

term  makes  negligible  contribution 
to  this  small  strain  constant  stress  result.  It  does,  however 
have  a  significant  role  when  S  /  0 ,  as  in  the  constant 
strain- rate  response  (35). 


Now  consider  an  initial  constant  shear  load  o_,  =  x 

21  o 


applied  over  (0,t^).  In  the  small  strain  approximation 
the  lead  order  invariants  and  stress  deviator  are  given  by 


(62)  and # (64),  and  neglecting  |y,e.,e2|  compared  to  unity  and 
noting  (e1,e2)  ^  0(y), 


I 

’  , 

I 

•  ,  .1 

r  .) 

|o  ^  Y 

l-2e1  y 

1 

—  li  1  ■», 

1  2  »  ( 

II 

,  p  = 

#  w  =  i  . 

, 

Y  l~2e2 

i 

2  ^ 

1  *•  1 

o 

-tr* 

1 

170) 


?  - 1 


4ytq  0 


0  -2yx 


,  tr(KI); 


Then 


/the  three  independent  components  of  (52)  are  to  lead  order: 


♦i*i  *  3*v  —  <lrov  +  2-  i°i  >  ■ 


*1*2  *  f*W 


=  y 
3  o 


2“  2  c  i 


(71) 


♦  jY  +  2ino  (e1+e2-jyr )  =  2(:o  -  u^y). 


Now  |ipTQ|  <<c  |<P^y|  for  the  model  (51), 

so  the  ip  term  in  the  shear  component  (71)  is  negligible, 

3 

but  if  | e^ , e 2 !  <<  Y  ,  the  v  and  4>  ^  terms  of  (71)^  2 

may  be  comparable.  By  (51) j ,  (71),,  |m^y|  <  0(xQ)  and 


n 


$^Y  =  0(tq)  ,  so  by  ( 51 )  ^  /  (71)  ^  we  see  that  4*^6^,  4>  1e^ 

_  2  _  _ 

are  given  by  terms  of  order  i~/K,  t  y,  t  e/y,  anc*  hence 

|(e^,e2)/y|  <■<  1  if  e  y  .  Initially,  the  jump  ratios 

[e^D/Cy]  i  [e2]/[Y3  are  negligible  because  of  the 

approximations  (63),  which  suggests  that  l®i,®2l  <<:  ^  °n 

(0,t^).  Figure  8(b),  (c)  show  tin.-  strain  response  for 

—  - 5—2 

a  =  0  and  a  =  0*5  when  \  Q  -  3  2  *  10  Nm  (same  J2  as 
axial  case) ,  determined  by  integrating  the  three  equations 
(71)  simultaneously;  indeed  and  e2  remain  negligible 

on  (0, t^)  . 

When  the  shear  load  is  removed  at  t  =  t^  there  are 
negligible  elastic  strain  jumps  followed  by  constant  strain 
on  (t-^,t2)  and  negligible  strain  jumps  on  reloading  at 
t  =  t2»  But  since  e^  and  e„,  remain  negligible  for 
0  £  t  ^  t2+  ,  subsequent  dependence  of  K^-3  is  symmetric 
in  the  axial  strain  changes  e^  e^(t2+)  and 
e2  -  e2(t2  +  ),  where  e1(t2+)  <?.,  (t2  +  )  =  0,  and  also  the 

respective  axial  loading  equations.  Hence  the  response  e^(t) 
to  constant  stress  °n  “  "  '’0  1  n  t  >  t2  is  the  same  as 

e2(t)  for  ° 22  ~  ~  °o  ;  that  is,  the  unloaded  configuration 
after  a  shear  cycle  is  symmetric  to  subsequent  axial  loading. 

Consider  the  response  to  =  -  oq  in  t  >  t2«  From  (62) 

2  ■> 

with  e2  =  -  ^e^,  K^-3  =  3e^  +  y "  ,  now  depending  on  y  ,  and 

by  (58) ,  B ^ 2  =  B2^  =  Y  ,  so  the  axial  equation  (65)  appliesr 
coupled  to  y  through  K^-3  ,  and  the  shear  component  of 
(52)  gives 


$3_Y  =  "  2w1y 


(72) 


i; 


to  be  integrated  siniu  1 '  ..in,- .  y  i  «  ::i  starting  values 
el  (t2  +  )  =  Y(t2  +  ).  When  ,  (■»  O)  ,  y  remains 

constant,  but  for  u  '■  »>  *.!•••  i  .  i  if  ieunr  shear  strain 

decay  during  the  axial  re  l  nd.:,:,  •  et  l.-ctinq  the  dependence 
on  the  strain  tensor.  Figui.  u  these  responses;  for 

a  =  0  and  a  =*  0*  5 .  him:  : .  •  » I  .•••  .its  for  >  ,  e2  are 
obtained  if  the  alternative  :  ■  ■  .iing  ■  ,)r)  -  a  is  applied. 

Note,  however,  if  alternat i ve  .tr  reloadings  =  7q 

and  -  t 0  are  applied,  the  response  will  be  different 

because  of  the  roles  of  the  cor i  -'spending  shear  strain  y^ 
and  y^  in  the  invariant. 

Concluding  remarks 

It  has  been  shown  that  a  differential  operator  relation 
defining  an  incompressible  viscoelastic  solid  can  describe  the 
observed  time  response  of  ice  to  applied  constant  load  and  to 
applied  constant  displacement  rate  in  uni-axial  stress. 
Furthermore,  both  strain-rate  and  stress-rate  terms  are 
necessary,  in  addition  to  strain  and  stress  terms,  to  reproduce 
both  these  responses.  By  comparison  there  is  no  strain 
dependence  in  the  analogous  incompressible  viscoelastic  fluid 
relation,  but  a  strain-acceleration  term  is  required,  so 
that  both  forms  of  relation  contain  three  independent  response 
coefficients  in  uni-axial  stress.  In  the  solid  model  it  is 
assumed  that  the  response  coefficients  depend  only  on  stress 
and  strain  (not  on  rates) ,  and  the  above  two  responses 
determine  only  two  combinations  of  the  three  response  coefficiei 
Thus  the  general  uni-axial  response  described  by  the  model  is 
not  fully  determined  by  these  two  responses,  in  contrast  to 
the  analogous  overdetermined  viscoelastic  fluid  model.  Further 
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investigation  of  the  construction  of  a  solid  model  is 
required.  The  solid  model  does,  however,  in  contrast  to 
the  fluid  model,  exhibit  strain  jumps  when  stress  jumps 
are  applied.  We  have  derived  the  general  jump  relations 
from  the  differential  tensor  relation,  and  directly  for  uni¬ 
axial  stress  from  the  reduced  uni-axial  stress  relation,  by 
assuming  that  they  are  the  limit  of  continuous  changes 
taking  place  in  a  decreasing  time  interval. 

Uni-axial  response  does  not,  of  course,  determine  the 
tensor  shape  of  the  differential  relation,  nor  the  dependence 
of  response  coefficients  on  stress  and  strain  tensors.  Test 
responses  involving  two  independent  stress  and  two  independent 
strain  components  are  necessary  to  determine  such  shape.  For 
illustration  we  adopt  simplifying  shape  assumptions  and  construct 
a  reduced  tensor  relation  compatible  with  idealised  uni-axial 
stress  response.  The  same  uni-axial  response  would  give  a 
different  tensor  relation  if  different  shape  assumptions  are 
made.  It  is  shown  in  general  that  the  response  from  a  pre- 
loaded  distorted  state  is  anisotropic,  even  when  the  initial 
configuration  is  isotropic.  Such  induced  anisotropy  is 
demonstrated  by  the  responses  of  the  reduced  model  to  loadingr 
unloading-reloading  cycles  when  the  reloading  is  applied  in 
different  directions.  Both  uni-axial  and  shear  loading  is 
considered  for  the  initial  load  stage.  Furthermore,  reloading 
which  repeats  the  initial  loading  shows  the  change  of  response 
which  is  caused  by  strain  dependence  in  a  solid  model,  not 
present  in  aifluid  model.  Correlation  of  test  data  with  solid 

I 

models  must  therefore  indentify  the  sample  configuration  with 
respect  to  an  undistorted  state.  If  ice  completely  relaxes  on 
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unloading,  given  sufficient  time,  then  a  test  sample  would 
have  the  original  configuration  on  formation,  isotropic  or 
anisotropic.  The  model  adopted  tor  illustration  exhibited 
relaxation  on  unloading. 
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Typical  response  for  constant  load  test  on  ice: 

a)  creep  curve, 

b)  strain-rate  v.  time. 

Typical  stress-strain  response  for  constant  displacement- 
rate  test  on  ice. 

Strain  relaxation  after  complete  unloadin-j: 

(i)  no  viscoelastic  relaxa: ion, 

(ii)  partial  viscoelastic  re  !  axe.  t  ion , 

(iii)  complete  viscoelastic  relaxation. 

Idealised  response  for  constant  Load  tests: 
a)  normalised  creep-curve. 


b)  -  f  (e)  —  f  (t)  . 


Calculated  stress-strain  rosir  use  for  four  constant 
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Response  to  load-unload-repeat  load  in  uni-axial  stress, 
all  a  . 

Response  to  load-unload-new  direction  reload  in  uni-axial 
stress,  for  a  =  0,  (b)  ,  :  =  o*r>  ic)  .  Comparison  of  e ^ 

and  in  t  >  reflects  the  induced  anisotropy  in 

states  at  t  =  t 

Response  to  shear  load-unload  followed  by  uni-axial  load 


for  a 


0 .  ( b )  ,  a  =  0  •  5  ( c )  . 


Figure  1.  Typical  response  Lm  constant  loud  test  on  ice 
(a)  creep  curve 
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(b)  strain-rate  v.  time. 


different  E(e): 
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Figure  7.  Response  to  load-unload-new  direction  reload  in  uni-axial 
stress,  for  a  =  0,  (b)  ,  a  --  0*5  (c)  .  Comparison  of  e2 
and  e,  in  t  >  t0  reflects  the  induced  anisotropy  in 
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Abstract 

The  finite  linear  viscoelastic  solid  and  single  integral 
representation  with  non-linear  dependence  on  history  are 
investigated  in  uni-axial  stress.  Both  integral  kernels  in 
the  stress  formulation  are  determined  by  single-step  constant 
strain  tests,  and  both  kernels  in  the  strain  formulation  are 
determined  by  single-step  constant  stress  tests.  Single  integral 
stress  and  strain  formulations  are  not  equivalent.  The  stress 
histories  required  to  maintain  constant  strain-rate  for  both 
models  are  determined  from  the  Vol terra  integral  equations  given 
by  the  strain  formulations  once  their  kernels  are  determined  by 
constant  stress  tests.  However,  known  constant  strain-rate 
response  does  not  determine  the  kernels.  Examples  are  presented 
to  show  that  variation  of  the  kernel  within  a  given  qualitative 
shape  can  lead  to  different  shapes  of  constant  strain-rate 
response,  so  that  both  constant  stress  and  constant  strain-rate 
tests  may  be  necessary  to  deduce  the  optimum  single  integral 
approximation,  in  preference  to  multi-step  stress  tests.  It  is 
shown  that  the  apparently  simpler  finite  linear  viscoelastic  model 
requires  a  far  lengthier  numerical  algorithm  to  solve  the  Volterra 
equation,  and  leads  to  non-unique  and  physically  unacceptable 
response,  in  comparison  with  the  more  flexible  non-linear  history 
dependence  which  yields  unique  acceptable  responses. 
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where  the  stress  deviator  S  is  defined  by 


S  =  o  +  pi  ,  p  =  -  |  tr  o  ,  (1.5) 

and  the  mean  pressure  p  is  a  workless  constraint  not  determined 
by  the  deformation.  Lockett  (1972)  suggests  that  a  given  pure 
shear  strain  history  (of  an  incompressible  material)  may  require 
isotropic  pressures  p,  in  contrast  to  the  shear  relations  (1.4) 
based  on  the  workless  constrain  postulate.  In  this  case,  the 
incompressibility  condition  (1.3)  should  be  adopted  as  a  good 
kinematic  approximation  (small  dilation  compared  to  typical 
shear  strains),  and  the  full  stress  relation  (1.1)  or  (1.2) 
applied . 

Theoretical  modelling  of  non-linear  viscoelastic  behaviour 
has  focussed  mainly  on  multiple  integral  expansions  of  the 
stress  functional  ,  (1.1),  with  kernels  depending  only  on  i 

time  lapse  t  -  t  (Green  and  Rivlin  1957,  1960;  Green,  Rivlin  ; 
and  Spencer  1959,  Pipkin  1964),  and  specific  representations 
for  isotropic  response  have  been  formulated.  Such  models  are 
usually  truncated  after  the  triple  integral,  but  still  require 
a  large  programme  of  tests  to  determine  the  kernel  functions. 
Lockett  (1965,  1972)  has  presented  theoretical  programmes  of 
multiple  step  constant  stress  creep  tests  sufficient  to  determine 
the  kernels  of  the  third  order  one-dimensional  and  three-dimensioni 
models,  with  corresponding  constant  strain  tests  for  a  third 
order  expansion  of  the  strain  functional  **(1.2).  The  third  fl 
order  expansion  of  the  strain  formulation  (1.4) 2  for  the 
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incompressible  case  is  analysed  further  by  Lockett  and 
Stafford  (1969) .  Motivation  for  the  truncated  stress  integral 
expansion  is  that  the  current  strain  and  weighted  integration 
of  past  strain  (fading  memory)  are  small.  Dill  (1975),  though 
the  expansion  is  an  exact  theory  which  may  describe  the  large 
strain  response  of  some  material.  Truncation  of  the  strain 
expansion  has  less  motivation,  and  truncations  of  both  expansions 
to  the  same  order  are  not  equivalent  formulations,  and  require 
different  numbers  of  kernels  in  the  incompressible  case 
(Lockett  1972) . 

An  alternative  expansion  in  a  small  strain  difference 
history  from  the  current  strain,  E(t)  -  E(t)  ,  leads  to 
integral  series  with  kernels  depending  on  the  fixed  strain 
point  E ( t) ,  (Dill  1975).  The  first  term,  linear  in  the 
strain/history ,  defines  a  finite  linear  viscoelastic  material 
(Coleman  and  Noll  1961)  ,  given  by 


a  (t)  =  K[ E  ( t)  ,  0]  E  ( t ) 


K[E  (t)  , 


J 

—  00 


t-T )E (x ) dT  , 


(1.6) 


where  K(E,t)  is  a  fourth  order  tensor  function,  zero  on 
t  <  0,  and  K  denotes  the  derivative  of  K(E,t)  with  respect 
to  its  second  argument,  supposed  continuous  in  both  arguments. 

The  first  term  of  a  similar  strain  expansion  in  stress  difference 
history  gives 


A  A 

E(t)  =  Cto (t) ,0]o ( t) 


C[o (t) ,t-x  ]o  (t)  dt  , 


(1.7) 


1 

j 


where  C[a,t]  is  a  fourth  order  tensor,  zero  on  t  <  0, 
and  C  denotes  derivative  with  respect  to  the  second  argument. 
It  is  convenient  to  describe  (1.7)  also  as  a  finite  linear 
viscoelastic  material. 

The  linear  viscoelastic  relations  for  infinitesimal  strain 
E  and  infinitesimal  rotation  R  -  1  are  given  by  (1.6)  and 

A 

(1.7)  with  o=o  and 

K(E,t)  =  K ( t)  ,  C (o , t)  =  C(t)  .  (1.8) 

In  this  case  (1.6)  and  (1.7)  are  equivalent  with  the  relaxation 
function  K ( t )  and  creep  function  C(t)  related  by  the  linear 
Volterra  integral  equation 

K (0)  C(t)  +  |  K’(t-T)  C(x)dT  =  1  H ( t )  (1.9) 

where  1  is  the  unit  fourth  order  tensor  and  H(t)  is  the 
Heaviside  step  function.  K(t)EQ  is  the  stress  response  to 
constant  strain  EQH(t),  and  C(t)oQ  is  the  strain  response 
to  constant  stress  oQH(t). 

From  the  finite  relations  (1.6)  and  (1.7),  K(EQ,t)E0  is 

the  stress  response  to  constant  strain  EQH(t),  and 
/\  ^  * 

C(aQ,t)ao  is  the  strain  response  to  constant  stress  oQH(t), 

but  there  is  no  integral  equation  analogous  to  (1.9)  which 

A 

determines  a  function  of  two  variables  C(o,t)  given  K(E,t), 
or  vice-versa.  Let  K(E  ,t)EQ  =  K(EQ,t)  be  given,  then  (1.7) 
gives 

■t 

Eq  =  CfMEo/f) /OlKtfo/t)  +  C[K(EQ,t)  ,t-T]K(E0,x)dT  (1.10) 

o 

which,  for  each  fixed  EQ,  is  an  integral  equation  to  determine 
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C[K(EQ,t) ,t]  for  0  <  t  <  t  ,  ts0.  As  different  Eq 
are  chosen,  the  stress  argument  K(EQ,t)  will  be  repeated 
at  different  values  t  ,  but  the  corresponding  solutions  of 
(1.10)  at  the  same  time  t  <  t  ^  will  not,  in  general,  be 
identical;  that  is,  there  is  no  C(o,t)  which  satisfies 
(1.10)  for  all  Eq.  This  becomes  more  evident  when  a  numerical 
algorithm  to  solve  the  one-dimensional  form  of  (1.10)  is  attempted 
Hence,  unlike  the  fully  linear  theory  (1.8),  (1.6)  and  (1.7) 

are  not  generally  equivalent  formulations.  Given  the  kernel 
C,  of  course,  (1.7)  is  an  integral  equation  for  the  stress 

A 

history  o(t)  necessary  to  support  any  prescribed  strain 
history  E(t)  ,  and  a  one-dimensional  numerical  algorithm  is 
presented  later.  Corresponding  results  follow  from  the  stress 
formulation  (1.6). 

Both  forms  of  multiple  integral  expansion  discussed  above  are 
motivated  by  a  strain  history  or  difference  history  remaining 
small,  with  corresponding  strain  expansions  less  firmly  based. 
Pipkin  and  Rogers  (1968)  point  out  that  strong  non-linearity  may 
then  require  a  high-order  multiple  integral  expansion,  with  a 
prohibitive  experimental  programme,  and,  furthermore,  very 
difficult  stress  and  deformation  analysis.  They  propose  that  the 
first  single  integral  should  account,  if  possible,  for  the  strong 
non-linearity  by  allowing  non-linear  dependence  of  the  kernel 
on  the  strain  or  stress  history,  with  successive  multiple 
integral  terms  viewed  as  refinements.  Both  stress  and  strain 
expansions  are  constructed  so  that  the  single  integral  describes 
exactly  the  response  to  single-step  constant  strain  or  constant 


stress  tests,  and  the  multiple  integral  terms  vanish  identically 
for  such  histories.  The  successive  multiple  integral  terms 
are  constructed  so  that  the  nth  order  kernels  are  determined 
by  an  n-step  constant  strain  or  constant  stress  test  with 
higher-order  terms  vanishing  identically.  Thus,  extra  terms 
can  be  added  to  refine  the  correlation  with  extra  multi-step 
response  without  disturbing  the  integral  series  already 
constructed,  in  strong  contrast  to  the  earlier  multiple  integral 
expansions  in  which  all  the  kernels  change  as  the  level  of 
truncation  changes.  The  single  integral  representations  are 


0(t) 


GfE(t),03  + 


t 

GfE(T)  , 


t-T Idr 


(1.11) 


E(t) 


/N 

J[  o (t) ,03  + 


t 


J 


Jl'O  (T  )  , 


t-T  3dt 


(1.12) 


where  the  relaxation  function  G[E,tl  and  creep  function 

A 

J[a,t]  are  second  order  tensors  depending  on  the  tensors 
E  and  a  respectively,  as  well  as  on  time.  G(E,t)  and 
J(a,t)  are  defined  to  be  zero  for  E  =  o  and  a  =  O 
respectively,  then  G(EQ,t)  is  the  stress  response  to  constant 
strain  EQH(t),  and  J[o  ,t]  is  the  strain  response  to 

A 

constant  stress  oQH(t).  An  integral  equation  corresponding 
to  (1.10)  is  given  for  the  constant  strain  response,  and 
as  before  the  kernel  G  does  not  determine  the  kernel  J  , 
nor  does  J  determine  G  ,  so  (1.11)  and  (1.12)  are 
not  equivalent.  The  creep-relaxation  relations  constructed 
by  Findley,  Lai,  and  Onaran  (1976)  are  based  on  iterative 


approximations  for  G  and  J  quadratics  in  the  histories. 

Lockett  (1972)  argues  that  successive  multiple  integral 
terms  in  the  Pipkin  and  Rogers  (1968)  scheme  are  not  readily 
determined,  since  the  number  of  arguments  in  the  kernels  increases 
rapidly  because  of  the  dependence  on  both  time  and  strain  or 
stress  history,  but  that  this  representation  is  a  useful  basis 
for  simplifications  such  as  kernels  seperable  in  time  and  strain 
or  stress  history.  Pipkin  and  Rogers  (1968)  have  used  single-step 
constant  uni-axial  stress  data  for  a  PVC  plastic  to  determine 
the  one-dimensional  form  of  the  single  integral  representation 
(1.2),  which  is  then  shown  to  predict  good  agreement  for  a 
five-step  stress  history  test.  Single  integral  representations 
for  polypropylene  and  polyethylene  also  predict  good  agreement 
for  two-step  stress  histories,  1  jt  for  a  polyurethane  foam  do 
not  agree  well  with  a  four-step  stress  test.  They  also  argue  that 
multi-step  stress  tests  provide  a  more  severe  trial  of  the 
single  integral  representation  than  smooth  stress  history 
response.  Clearly  single  integral  representations  are  the  most 
tractable  for  stress  and  deformation  analysis  and  require  the 
least  extensive  experimental  programme,  and  the  forms  (1.11) 
and  (1.12)  which  incorporate  non-linear  history  dependence 
offer  most  flexibility,  with  some  good  response  correlations 
already  demonstrated  by  Pipkin  and  Rogers  (1968) . 

We  now  inve.'  gate  the  strain  representation  (1.12)  and 
the  corresponding  ..  ..ite  linear  representation  (1.7)  in  uni-axial 
stress.  Given  the  scalar  kernels  J  and  C  defining  the 
same  single-step  constant  stress  response,  the  respective  stress 
histories  required  to  maintain  constant  strain-rate  are  determined 
by  Volterra  integral  equations.  Distinct  numerical  algorithms 


are  necessary.  It  is  found  that  the  apparently  simpler  finite 
linear  viscoelastic  relation  requires  a  considerably  lengthier 
and  less  satisfactory  algorithm  than  the  general  non-linear 
relation,  which  follows  also  for  other  loadings.  Furthermore, 
in  examples  with  C  determined  by  a  smooth  constant  stress 
response,  the  constant  strain-rate  stress  history  for  the  finite 
linear  viscoelastic  relation  becomes  non-unique  and  all  branches 
become  physically  unacceptable.  In  contrast,  the  general 
non-linear  relation  yields  a  unique  physically  acceptable  solution, 
and  accuracy  and  stability  of  the  numerical  algorithm  is  tested 
by  application  to  the  linear  viscoelastic  relations  for  which 
exact  solutions  can  be  derived.  However,  variation  of  the 
kerenel  J  within  a  given  qualitative  shape  can  lead  to  widely 
differing  shapes  of  stress  history  for  constant  strain-rate. 

The  same  effect  is  also  exhibited  by  solutions  of  the  linear 
viscoelastic  relations  when  similar  kernel  variation  is 
investigated,  and  is  not  peculiar  to  the  chosen  non-linear 
kernel.  Thus,  when  constant  strain-rate  response  is  known  as 
well  as  constant  stress  response,  the  "best"  single  kernel 
is  not  necessarily  the  exact  single  step  constant  stress  response, 
as  in  the  Pipkin  and  Rogers  (1968)  scheme,  but  rather  that  required! 
to  approximate  both  responses.  That  is,  the  trial  of  a  single 
integral  representation  should  be  multi-type  tests  rather  than 
multi-step  tests  of  the  same  type.  The  constant  strain-rate 
response  is  a  smooth  stress  history,  which,  in  our  examples, 
does  provide  significant  information.  Our  inverse  methods  do 
not  show,  however,  how  the  constant  strain-rate  response  can  be 
used  directly  to  "optimise"  the  kernel  J  ,  nor  to  construct 
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double-integral  terms  to  refine  the  model,  in  contrast  to 
the  Pipkin  and  Rogers  (1968)  scheme,  but  have  revealed  the 
possible  failure  of  models  constructed  by  single-type  tests. 


The  non-linear  viscoelastic  relations  (1.11)  and  (1.12)  become 


a  (t) 

e  (t) 


GCe (t) ,0]  + 


GCe (t ) , t-i ]dx 


JCa(t) ,0] 


+ 


■t 

j[cr  ( r  )  ,  t-T  ]dx 

* 

—  oo 


(2.5) 

(2.6) 


The  kernels  K(t),  C(t)  in  the  linear  viscoelastic  theory, 

depending  only  on  time,  are  determined  by  the  response  to  any 

applied  stress  or  strain  history.  However,  dependence  of  K  on 

e(t)  in  (2.2)  and  of  G  on  e(x)  in  (2.5),  implies  that  K 

and  G  can  be  determined  uniquely  only  by  families  of  constant 

strain  tests,  e(t)  =  eQH(t)  with  a  sequence  of  eQ  to  cover 

the  required  strain  range.  The  stress  response  is  then  precisely 

K[eo,t]eo  and  G(eQ,t)  respectively,  described  as  the 

relaxation  function.  Similarly,  the  creep  function, or  strain 

response  to  constant  stress  o(t)  Ho  H(t)  ,  is  C[o  ,t]o 

o  o  o 

and  J[cro,t]  from  (2.3)  and  (2.6)  respectively.  Other  smooth 
strain  or  stress  histories,  including  constant  strain-rate 
e(t)  =  rtH(t)  or  constant  stress-rate  o(t)  =  stH(t),  do  not 
lead  to  unique  kernels  unless,  of  course,  the  particular  single 
integral  representation  is  an  exact  description  of  the  class  of 
smooth  responses  chosen.  Furthermore,  there  is  no  simple  correlate 
of  the  function  K[e(t),t-T]  ,  t  <  t  ,  nor  of  the  function  j 

GCe (t), t-t]  ,  t  <  t  ,  with  the  stress  response  o(t)  for  a  1 

general  history  e(t),  and  corresponding  conclusions  hold  for  the  I 
strain  formulations  (2.3)  and  (2.6).  Also,  the  non-uniform  stress  I 


response  a(t)  to  constant  strain  does  not  determine  the 
kernel  C  in  (2.3)  nor  the  kernel  J  in  (2.6),  and  the 
non-uniform  strain  response  e(t)  to  constant  stress  does  not 
determine  K  in  (2.2)  nor  G  in  (2.5).  For  multiple  integral 
terms,  only  multi-step  constant  strain  or  multi-step  constant 
stress  tests  allow  direct  correlation  with  the  kernels  of  a 
stress  or  strain  formulation  respectively,  and  the  two  formulation* 
cannot  be  related. 

Now  suppose  that  the  kernels  C  and  J  have  been  1 

determined  by  families  of  constant  stress  tests,  so  that  the 
representations  (2.3)  and  (2.6)  describe  the  same  constant 
stress  response.  The  Volterra  integral  equations  (2.3)  and 
(2.6)  for  the  respective  stress  histories  a (t'  corresponding 
to  a  prescribed  strain  history  e(t)  can  be  solved  numerically 
to  investigate  the  different  responses  predicted  by  the  two 
representations.  We  will  present  the  algorithms  and 
illustrations  for  constant  strain-rate  eft)  =  rt  as  an  example 
of  a  smooth  loading  history  used  in  practice.  Given  K  and  G 
from  constant  strain  tests,  the  Volterra  integral  equations 
(2.2)  and  (2.5)  similarly  determine  the  responses  to  an  applied 
stress  history  o(t). 

First  consider  the  non-linear  relation  (2.6)  with 
e(t)  =  rtH(t)  ,  where  r  is  a  constant,  so  that  o(t)  =  0 


for  t  <  0 


Define 


where  6  is  a  small  time  interval,  and  interpret  t 
By  (2.6) 


0+ . 


o  =  J[oo,0]  .  (2.8) 

If  J[o,0]  >  0  for  all  a  >0;  that  is,  there  is  a  positive 
strain  jump  at  t  =  0  when  a  positive  stress  jump  is  applied 
which  is  the  typical  response  of  a  solid;  then 

=  0  (2.9) 

and  the  stress  history  starts  smoothly.  If,  however,  the  initial 
strain  jump  is  negligible  compared  to  typical  creep  strains, 
so  J[o,0]  =0  is  adopted  as  an  approximation,  then  oq  is 
not  determined  by  (2.8),  and  may  be  non-zero.  Applying  the 
trapezoidal  rule  over  the  first  time  interval,  (2.6)  gives 

r<5  =  JCa1,0]  +  -^6  { J[  oq  ,  6  ]  +  3(^,0]}  .  (2.10) 

Now  if  J[a,0]  =  0  ,  3[a,0]  >  0  ,  the  limit  6  -*•  0  gives 

j[oQ,0]  =  r  .  (2.11) 

We  assume  that  the  initial  creep  rate  3[a,0j  at  constant 
stress  o  increases  with  a  ,  so  (2.11)  yields  a  unique 
solution  aQ  for  each  r  .  Now  oQ  is  determined  by  (2.9) 
or  (2.11),  and  in  both  cases  (2.10)  is  an  implicit  equation  for 
o^.  Continued  application  of  the  trapezoidal  rule  over  successiv 
time  intervals  gives  for  n  =  2,3,... 
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nr 5  =  J[on,0]  +  j  f>  J[c  ,0] 


x  •  n-i  . 

+  <S{±  J [a  ,  t  ]  +  Z  J[a. ,t  ]}  , 

2  on  l  n-i 


(2.12) 


which  is  an  implicit  equation  for  c>n  once  c ^  (i  =  0,1,... n-1 
are  determined  by  previous  steps,  and  occurs  in  only 

two  terms . 

The  algorithm  (2.10),  (2.12),  with  (2.9)  or  (2.11),  can  be 
applied  to  the  linear  viscoelastic  relations  by  setting 
J[cr(T),t]  =  J(t)  o  (t)  in  (2.6)  so  that  JCa^t^]  becomes 
J(tj)°i.  An  exact  solution  of  (2.6)  can  then  be  derived 
(by  Laplace  transforms)  for  exponential  type  creep  functions, 
and  examples  have  shown  the  time  interval  magnitude,  relative 
to  the  time  scale  of  the  creep  function,  necessary  for  the 
algorithm  to  yield  an  accurate  solution.  Further,  with  the 
assumption  that  J[o,0]  increases  with  a  ,  and  the  assumption 
J[o,0]  increases  with  a  ,  describing  expected  response, 

(2.10),  (2.12)  yield  unique  solutions  ai'a2''"'  * 

For  the  finite  linear  viscoelastic  relation  (2.3), 
constant  strain-rate  given  by  e(t)  =  rtH(t)  yields  the  initial 
condition 

0  =  C[oo,0]oo  ,  (2.13) 

with  solution  oq  =  0  if  CTo,0]  >  O  for  all  o  ,  and  the 
first  step  relation 

r<5  =  C[o.  ,0]o.  +  ^{CCo,,  6]o  +  Cfo.fOla.}.  (2.14) 


When  C[o,0]  =  0 


C[o,0]  >  O,  the  limit 


6  O  gives 


C[oo,0]oo  =  r 


(2.15) 


As  before  we  assume  C[a,0]  increases  with  a  ,  then  (2.15) 
yields  a  unique  solution  oQ  for  each  r  .  Again,  (2.14) 
is  an  implicit  equation  for  ,  but  here  all  three  terms  in 
C  have  argument  ,  and  in  particular  the  term  C[o^,6] 
involves  the  kernel  evaluated  at  the  unknown  and  time 

6^0.  Over  successive  time  intervals  the  algorithm  gives 
for  n  =  2,3,... 


nr6  =  C[ a  ,0 ]<j  +  (  Cfo  ,t  ]c>  +  CC a  ,0]o  } 

n  n  n  n  o  n  n 

n-1  . 

+  6  I  CTo  ,t  .]o, 

i=l  n  n-i  i 


(2.16) 


which  is  an  implicit  equation  for  once  cr^  (i  =  0,1,...  n-1) 

are  determined  by  previous  steps.  However,  every  term  in  C 

involves  the  unknown  argument  and  n  terms  involve 

CCo  ,t  .]  evaluated  at  t,,t,...  t  i  0  .  Thus,  solution 
n  n— i  il  n 

of  the  implicit  equation  (2.16)  at  each  step  becomes 

increasingly  lengthy  as  n  increases,  in  contrast  to  (2.12) 

where  o  occurs  in  only  two  terms,  and  those  evaluated 
n 

only  at  t  =  O.  Properties  of  C[a,0]  no  longer  allow  any 

conclusion  about  each  solution  o  since  evaluations  at 

n  n 

and  each  time  step  occur.  The  collapse  of  (2.16)  for  the  linear 
viscoelastic  relation  is  too  dramatic  to  provide  a  sensible 
test  of  the  numerical  algorithm,  but  solutions  have  been  j 

repeated  with  different  time  intervals  to  check  consistency.  I 
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Our  illustrations  show  that  non-uniqueness  can  indeed  occur 
for  smooth,  physically  sensible,  creep  functions  C,  so  not 
only  is  (2.16)  unattractive  for  numerical  methods,  but  the  finite 
linear  viscoelastic  parent  relation  (2.3)  can  imply  non-sensible 
response.  This  situation  must  also  arise  for  less  simple 
loadings . 

3.  Illustrations 

Constant  stress  and  constant  strain-rate  in  uni-axial 

compressive  stress  are  the  common  tests  used  to  determine  the 

mechanical  properties  of  ice  at  constant  temperature,  reviewed 

by  Mellor  (1980) .  Consider  a  and  e  to  be  positive  in 

compression.  While  detailed  results  over  the  large  temperature 

range  of  interest  (230K  -*■  273K)  have  not  been  determined, 

the  main  features  for  strains  up  to  about  0*05  are  known,  and 

typical  responses  are  shown  in  Fig.  1.  At  constant  stress  a  , 

Fig.  la,  there  is  an  initial  small  elastic  strain  jump 

e  =  o/E  ,  where  E  is  the  Young’s  modulus  at  zero  stress, 
e  o  o 

10  -2 

of  order  10  Nm  ,  followed  by  a  primary  decelerating  creep 

(e  <  0)  until  time  tm(o)  when  e  =  0  ,  e  =  rm(o), 

e  =  em(o)  »  then  a  tertiary  accelerating  creep  (e  >0).  The 

minimum  strain-rate  r  (a),  often  referred  to  as  secondary 

or  steady  state  creep,  is  significantly  non-linear  in  a  , 

and  highly  temperature  dependent.  A  close  least  squares  fit  to 

laboratory  data  at  273.13K  over  the  stress  range 
5  -2 

0  s  as  9  *10  Nm  derived  by  Smith  and  Morland  (1981)  is 
r  <o)  =  o(0*3336  +  0*  320c ?  +  0*0296o“) 


(3.1) 
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i 


5  -  2 

where  a  is  measured  in  units  10  Nm  and  strain-rate  in 
units  a  The  strain  at  minimum  strain-rate,  em(o),  is 
approximately  0*01,  independent  of  a  over  this  stress  range, 
and  much  larger  than  the  elastic  strain  e  (o).  Figure  lb 
shows  the  corresponding  strain-rate  at  constant  stress.  The 
typical  response  to  constant  strain-rate  e  =  r  is  shown 
in  Fig.  lc,  starting  smoothly  at  zero  stress  with  a  increasing 
to  a  maximum  a„(r )  at  time  t„.(r),  then  decreasing.  The 
strain  at  tM(r),  eM(r)  =  rtM(r)  '  is  a-*-so  approximately  0*01 
for  a  moderate  range  of  r  .  While  the  creep  at  constant 
stress.  Fig.  la,  is  monotonic,  the  decelerating  and  accelerating 
stages  shown  more  clearly  in  Fig.  lb  are  not  common  in  other 
rheological  models.  There  is  a  corresponding  non-monotonic 
stress  response  at  constant  strain-rate,  Fig.  lc,  but  we  find 
that  this  does  not  follow  from  (2.6)  for  all  creep  functions 
with  the  above  properties . 

An  idealised  creep  model  for  the  respresentation  (2.6) 
which  reflects  the  qualitative  features  of  ice  deformation, 
adopted  by  Morland  and  Spring  (1981)  and  Spring  and  Morland  (1982) 
to  illustrate  differential  operator  relations,  is  given  by 


j[0,t]  =  rm  (a )  { A  -  Be~bt/tfn  +  De  dt/tm)  , 


(3.2) 


J[o,t] 


>t 

j[o,t]dt  +  0/Eo  , 
o 


(3.3) 


with  r  (o)  defined  by  (3.1)  and  E  = 
m  o 

making  the  approximation  em(o)  2  0*01, 


9  x  lO^Nm  Then 
(3.2)  and  (3.3)  give 
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°*01  "  °/E0  ,  B ,  -b  ,  %  D,  -d  ,  ,  ,, 

tm(0)  =  ~ mt) -  '  a  =  A  +  b(e  -U-gte  -1),  (3.4) 

m 


where  the  constant  a  must  be  positive.  Also 


J[a,0]  =  a/Eo  >  0,  j[a ,0]  =  r m(o){A  -  B  +  d}  ,  (3.5) 


so  J[o,0]  and  j[o,0]  increase  with  a  provided  that  A  -  B  +  D  >  0. 

The  algorithm  (2.12)  has  been  used  to  determine  the  stress  history 
for  three  constant  strain-rates  r  =  0*5a  \  la  \  2a  \  and 
the  three  parameter  sets  for  (3.2)  and  (3.4)  listed  in  Table  1. 


A 

B 

D 

b 

d 


I  II 

2*5  2* 

3*6664  1* 

10*7206  15* 

0*6706  0* 

5b 


III 

5 

2*5 

5660 

1*5164 

1685 

29*2772 

0380 

0*0099 

200b 

lOOOb 

Table  1.  Parameters  for  idealised  creep  functions 


These  all  give  the  same  ratio  t^/t^  =  1*7  ,  where  t^.  is  the  time 
to  the  inflexion  point  (J  =  0)  in  the  tertiary  creep,  and  satisfy 
a  >  0  and  J[o,0],  J[o,0]  increasing  in  a  .  From  (3.2)  -  (3.4), 


0*01  -  o/E, 


J(a,t)  =  ~ 
Eo 


-{A 


_t_ 

t 


m 


+  i  (e 


-bt/t, 


m-l)  -  f(e 


-dt/t, 


m 


-1)}  (3.6) 


is  a  function  of  t/tm(o)  alone  when  o/E0  is  neglected.  Figure 

2  shows  e  and  e/r(a)  as  functions  of  t/t  (o)  for  the 

m  m 

parameter  sets  I,  II,  III;  each  has  the  required  qualitative 
shape . 


The  distinction  between  sets  I,  II,  and  III,  is  seen  in 
the  stress  histories  they  predict  at  constant  strain-rate,  shown 
in  Fig.  3,  as  functions  of  e  =  rt.  With  set  I  there  are  large 
repeated  oscillations,  with  set  II  there  is  a  small  stress  rise 
after  relaxation  before  the  steady  response,  and  with  set  III 
the  required  stable  shape  corresponding  to  Fig.  lc  is  obtained. 
If  a  linear  viscoelastic  relation  with  related  creep  function 


Bt 


J(t)  =  ~~  +  j { At  + 
o 


m 


(e 


-bt/t, 


m 


Dt 


-  1  )  - 


m 


(e 


-dt/t, 


m 


-  1)  }  , 


(3.7) 


where  j  and  t  are  constants,  is  solved  exactly  for  the  constant 
m 

strain-rate  response,  then  oscillations  during  the  relaxation 
phase  arise  if 

X  .  [ A (b+d)  -  Bd  +  Db]-  <  J  _  ,3.8, 

4 Ad  (A  -  B  +  D) 


when  the  elastic  response  1/Eq  is  neglected.  For  the  parameter 
sets  I,  II,  III,  X  =  0*11,  1*3,  and  3*4  respectively,  so  that 
small  oscillations  still  occur  during  relaxation  in  the  non-linear 
model  at  X  =  1*3,  but  not  at  X  =  3«4.  The  criterion  (3.8) 
provides  a  rough  guide  to  the  likelihood  of  relaxation  oscillations, 
These  illustrations  show,  however,  that  creep  functions  of  the 
correct  qualitative  shape  can  lead  to  distinct  shapes  of  constant 
strain-rate  response  from  the  representation  (2.6),  so  that  an 
optimum  approximation  for  the  kernel  J  must  take  account  of 
constant  strain-rate  response.  The  above  procedure,  of  course, 
provides  only  an  inverse  approach,  and  a  direct  method  has  yet 


to  be  constructed. 
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The  same  constant  stress  response  (3.2)  -  (3.4)  is  described 
by  the  finite  linear  viscoelastic  relation  (2.3)  with 

C[a,t]  =  J[c,t]/a,  C[  o ,  t]  =  j[o,t]/o  ,  (3.9) 

which  are  bounded  as  a  -*■  0  .  Application  of  the  lengthy 
algorithm  (2.16)  with  the  parameter  set  II  and  r  =  la  ^  leads 
to  a  non-unique  stress  long  before  the  peak  stress  given  by  the 
relation  (2.6)  is  reached,  and  the  time  at  which  this  is  observed 
numerically  depends  on  the  time  interval  6  .  Table  2  provides 
an  outline  of  the  numerical  results  obtained  with  time  intervals 
6=1  hour,  2  hours,  10  hours,  corresponding  to  strain  increments 
1*142  x  10  2*284  x  lo  ^ ,  1*142  x  lo  ^  respectively  at 

r  =  la  \  and  the  results  from  the  non-linear  relation  algorithm 
(2.12)  with  6=2  hours  for  comparison.  Once  non-uniqueness 
occurs,  there  are  three  roots  of  the  implicit  equation,  and  each 
branch  yields  three  roots  at  the  next  step,  so  no  physically 
sensible  solution  can  be  deduced.  The  Table  shows  only  the 
continuation  of  th°  minimum  and  maximum  root  branches  of  the 
first  non-unique  solution.  Until  non-uniqueness  occurs  the 
solutions  for  different  6  are  consistent,  and  close  to  the  general 
non-linear  solution  as  expected  for  such  a  short  initial  time 
range.  This  dramatic  non-uniqueness  property  of  the  finite  linear 
viscoelastic  integral  equation  must  extend  to  other  prescribed 
strain  histories,  and  the  representation  (2.3)  would  appear  to  have 
no  general  validity.  The  same  situation  must  arise  in  the  usual  j 
stress  formulation  (2.2)  with  some  smooth  kernels  K[e,t]  when 
the  integral  equation  for  constant  stress-rate  response  is  solved. 


hours 


a  C 10  Nm  15=1  hour 


6=2  hours 


6=10  hours 


6=2  hours 

(2.12) 


0 

0*181 

0*181 

0*181  ; 

1 

0*183  ; 

1 

1 

2 

0*186 

0*185 

i 

1 

3 

0*188 

4 

0*191 

0*191 

5 

0*194 

6 

0*197 

0*197 

7 

0*200 

8 

0*203 

0*203 

9 

0*207 

10 

0*210 

0*210 

0*206 

11 

0*214 

! 

12 

-  0*018  0*218  0*054 

0*218 

13 

0*068  0*222 

14 

-  0*200  0*226 

0*266 

15 

0*808  0*254 

16 

-  0*351  -  0*093 

-0*021  0*236 

17 

1*127  1*183 

18 

-  0*630  -  0*627 

0*876  0*241 

19 

-  0*037  -  0*033 

20 

1*511  1*496 

-  0*364  0*393 

0*251 

22 

I'  'i' 

1*111  -  0*169 

24 

o  o  1 

-  0*505  1*194 

min  max 

O  0 

26 

-  0*063  -  0*493 

min  max 

■  i 

28 

1*381  -  0*092 

i  i 

30 

-  0*415  0*012 

-0*069  0*322 

40 

f  f 

0*903  0*442 

50 

o'  o' 

-0*196  0*153 

min  max 

Table  2.  Algorithm  (2.16)  solution  for  different  time 
intervals  6  . 


•181 

•191 

•202 

0*213 

0*226 

0*236 

0*252 

0*267 

0*316 
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4 .  Concluding  remarks 

The  finite  linear  viscoelastic  representation  (2.3)  or  (2.2) 
appears  to  have  no  validity  for  significantly  non-linear  response 
with  the  features  described  above,  namely,  a  monotonic  strain 
response  to  constant  stress  with  an  inflexion  point,  or  analogous 
stress  response  to  constant  strain.  The  general  non-linear  single 
integral  representation  (2.6),  and  presumably  (2.5),  yields 
integral  equations  with  stable  numerical  algorithms  and  unique 
solutions  in  the  examples  treated.  However,  the  distinct  shapes 
of  constant  strain-rate  response  which  can  arise  for  kernels  of 
the  same  shape  determined  by  constant  stress  response  suggests 
that  multi-type  tests,  rather  than  multi-step  tests  of  the  same 
type,  are  necessary  to  obtain  a  valid  single  integral  approximation 
A  direct  procedure  for  correlating  the  kernel  with  multi-type 
response  has  yet  to  be  formulated,  and  may  prove  a  major  hurdle, 
but  the  inverse  approach  used  in  the  above  examples  would  offer 
a  trial  and  refinement  method. 

Acknowledgement 

This  research  has  been  supported  by  the  United  States  Army 
European  Research  Office  in  collaboration  with  the  Cold  Regions 
Research  and  Engineering  Laboratory,  Hanover,  New  Hampshire. 


23 


References 

Coleman,  B.  D.  and  Noll,  W.  1961.  Foundations  of  linear 
viscoelasticity,  Rev.  Mod.  Phys . ,  3j3,  239-249. 

Dill,  E.  H.  1975.  Simple  materials  with  fading  memory,  Chap.  4 
of  Continuum  Physics  II  (Ed.  Eringen,  A.C.),  283-403. 

Academic  Press,  New  York. 

Findley,  W.  N.,  Lai,  J.  S.  and  Onaran,  K.  1976.  Creep  and 
Relaxation  of  Nonlinear  Viscoelastic  Materials,  367pp, 
North-Holland  Pub.  Co.,  Amsterdam. 

Green,  A.  E.  and  Rivlin,  R.  S.  1957.  The  mechanics  of  nonlinear 

materials  with  memory.  Part  I,  Arch.  Rat.  Mech.  Anal.  1,  1-21. 

Green,  A.  E.  and  Rivlin,  R.  S.  1960.  The  mechanics  of  nonlinear 
materials  with  memory.  Part  III,  Arch.  Rat.  Mech.  Anal.,  4, 
387-404. 

Green,  A.  E. ,  Rivlin,  R.  S.  and  Spencer,  A.  J.  M.  1959.  The 

mechanics  of  nonlinear  materials  with  memory,  Part  II,  Arch. 
Rat.  Mech.  Anal.,  3,  82-90. 

Lockett,  F.  J.  1965.  Creep  and  stress  relaxation  experiments 
for  nonlinear  materials,  Inst.  J.  Eng.  Sci.,  3,  59-75. 

Lockett,  F.  J.  1972.  Nonlinear  Viscoelastic  Solids,  195pp, 
Academic  Press,  London. 

Lockett,  F.  J.  and  Stafford,  R.  1969.  On  special  constitutive 

relations  in  nonlinear  viscoelasticity,  Int.  J.  Eng.  Sci.,  7,  j 
917-930.  j 

Mellor,  M.  1980.  Mechanical  properties  of  polycrystalline  ice, 
in  Proc.  IUTAM  Symposium  :  Physics  and  Mechanics  of  Ice, 
Copenhagen  1979  (Ed.  Tryde,  P.),  117-128,  Springer-Verlag, 
Berlin. 

Morland,  L.  W.  and  Spring,  U.  1981.  Viscoelastic  fluid  relations 
for  the  deformation  of  ice.  Cold  Regions  Sci .Technol . ,  4,  j 
255-268.  j 

Pipkin,  A.  C.  1964.  Small  finite  deformations  of  viscoelastic 
solids,  Rev.  Mod.  Phys.,  36,  1034-1041. 

Pipkin,  A.  C.  and  Rogers,  T.  G.  1968.  A  nonlinear  integral 

representation  for  viscoelastic  behaviour,  J.  Mech.  Phys. 
Solids,  16,  59-72. 


i 


24 


Smith,  G.  D.  and  Morland,  L.  W.  1981.  Viscous  relations  for 
steady  creep  of  polycrystalline  ice.  Cold  Regions  Sci. 
Technol.,  5,  141-150. 

Spring,  U.  and  Morland,  L.  W.  1982.  Viscoelastic  *=•  1  id  relations 
for  the  deformations  of  ice.  Cold  Regions  Sci.  _;Chnol. 

(in  press)  . 


Figure  Captions 


1.  Typical  response  of  ice  in  uni-axial  compressive  stress 

(a)  strain  at  constant  stress, 

(b)  strain-rate  at  constant  stress, 

(c)  stress  at  constant  strain-rate. 

2.  Strain  and  strain-rate  at  constant  stress  for  the 

three  models  I  ( - ),  II  ( - ),  III  ( - ). 

Stress  at  constant  strain-rates  r  =  »5a  \  la-\  2a~^, 
for  the  models  I,  II,  and  III. 
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Strain  and  strain-rate  at  constant  stress  for 
the  three  models  I( - ),  II  ( - ),  in)-.-. 


Figure  3.  Stress  at  constant  strain-rates  r 
2a  ,  for  the  models  I,  II,  and  III 


